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Geometry of spectral curves and all order dispersive 

integrable system 
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O ■ Abstract 
(N 

CN ' We propose a definition for the tau function and a Baker-Akhiezer spinor kernel of an 

integrable system whose times parametrize slow deformations (at a speed of a given 

P_ii algebraic plane curve. This definition is a full asymptotic series in involving theta 

^ ! functions. The large limit of this construction reproduces the algebro-geometric solutions 

■ of the multi-KP equation. We check that our tau function satisfies Hirota equations to the 

first two orders, and we conjecture that they hold to all orders. The Hirota equations are 
equivalent to a self-replication property for the Baker-Akhiezer spinor kernel. We analyse 
I its consequences, namely the reconstruction of a isomonodromic problem given by a Lax 

^ [ pair, and the relation between "correlators" , the tau function and the Baker-Akhiezer spinor 

^ ! kernel. This construction highlights the bridges between symplectic invariants, integrable 

O ■ hierarchies and enumerative geometry. 



><j ■ 1 Introduction to integrability 



Integrable systems play a particularly important role in physics, because they are 
"exactly solvable", at least in principle. Under a suitable change of variables, they can 
be brought into the form of a linear constant motion with constant velocity in a multi- 
dimensional torus. However, there is a difference between exactly solvable and exactly 
solved, and the study of integrable systems is still a very active branch of theoretical 
physics and mathematical physics. See the books [DJMOO, BBT02] for an introduction. 

There exists strong links between geometry and integrability [Kri77, SS83, SW85, 
Hit87, ManSS, Dub92, Kri94, IKF95]. Besides, the interplay with algebraic geometry 
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is an important issue. Indeed, generating functions which encode the enumerative ge- 
ometry in various kinds of complex manifolds have been shown to satisfy integrable 
equations. It happens for instance for generating functions for intersection numbers 
[Wit91, Kon92, Giv03], Hurwitz numbers and Hodge integrals [OkoOO, Kaz09], or 
Gromov-Witten invariants in toric Calabi-Yau threefolds [ADK+06, Brill]. 

There are several definition of what is an integrable system. In this article, we 
speak of an integrable system whenever some Hirota equations [Hir71] are present. An 
important class of integrable systems arise from Lax systems. The locus of eigenvalues 
of the Lax operator is an algebraic curve called the spectral curve. Here, conversely, we 
shall deal with the ways to build an integrable system (and in fact a Lax pair) starting 
from an arbitrary algebraic spectral curve, the so-called reconstruction problem. This 
has been addressed, in a dispersionless limit, in several earlier works (especially [Kri77] 
and [BG07]), on which we will comment along this article. Here, we propose a way to 
add the dispersionful part, using the algebraic invariants introduced in [EO07a]. 

Let us recall first generalities about Lax pairs and integrable systems. 

1.1 Hamiltonian system and Liouville integrability 

Consider a classical Hamiltonian system, i.e. a set of coordinates qi{t) and their con- 
jugated momenta Pi(t), obeying the Hamilton equations of motions: 

— p - -— —q-— (1 1) 

dt * dqi ' dt ' dpi 

where H is the Hamiltonian. This can be written in vector notation: 

Q = (gi, . . . . . . ^ = {H,Q} (1.2) 

where {, } is the Poisson bracket such that {pi, qj} = Sij. 

The system is "Liouville integrable" if there exist n independent conserved quanti- 
ties Hi, H2, . . . , Hn (by convention Hi = H): 

{H,,H,} = (1.3) 

In this case, it is possible to find a symplectic change of coordinates in the phase space, 
called "action-angle" variables, in which the equations of motions reduce to linear 
constant motions in a torus (the Jacobian, i.e. Ui is a periodic variable): 

dui dvi 

IF = dF = ° ^'-^^ 

The main problem is to find explicitly the conserved quantities H^, the action variables 
Vi, and the angle variables Ui as functions of the initial system of coordinates qi,Pi. 

A way to achieve this, was found in terms of a Lax pair, which we review below, 
and all the quantities above can be expressed in terms of algebraic geometry. 



2 



1.2 Lax matrix formulation 

Imagine that one can group the coordinates g^, pj together, using formal series 
Lij{x, Q), i,j = 1, . . . , d of a formal variable x, which form a. d x d matrix L{x, Q). 
This means that the coefficients of the expansion of Lij{x,Q) in powers of x, are 
functions of the coordinates q,p. Imagine that there exists a matrix M{x,Q), such 
that the equations of motion Eqn. 1.2 can be rewritten as a matrix equation: 

^ L{x, Q{t)) = [L{x, Q{t)), M{x, Qm (1.5) 

This equation is called a Lax equation, L is called the Lax matrix, x is called the 
spectral parameter, and (L, M) is called a Lax pair. It is easy to see that the Lax 
equation implies that for any x, the eigenvalues of L{x, Q(t)) are conserved quantities, 
and by expansion into powers of x, we obtain all conserved quantities. 

It may not seem obvious that the Hamilton equations of motion Eqn. 1.2 of a 
Liouville integrable system can always be written in this way, but indeed a Lax pair'^ has 
been found for all classical integrable systems (see for example [BBT02] and references 
thereof). 

1.3 The spectral curve 

The Lax equation (Eqn. 1.5) implies that the eigenvalues of L{x,Q{t)) do not depend 
on the time t, they are conserved. In particular, this means that the characteristic 
polynomial S{x,y) is independent of time: 

S{x,y) = det{y id- L{x,Qm (1.6) 

The locus of zeroes of £{x,y), i.e. the eigenvalues of L, is called the "spectral curve" 
S: 

S = {{x,y)\S{x,y)=0} (1.7) 

By expanding S{x, y) into powers of x and y, one can find all the conserved quantities, 
and thus, knowing the spectral curve, one knows the conserved quantities. The non- 
conserved quantities, i.e. the angle variables, are contained in the eigenvectors of L. 
It was shown [Kri77, BBT02], that the eigenvectors of L, must be Baker- Akhiezer 
functions on the spectral curve. Before explaining what it means, we need to describe 
more geometry of the spectral curve, i.e. more algebraic geometry. 

^ Given a Hamiltonian system, there is neither existence, nor unicity, of an equivalent Lax system. 
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1.4 Outline of the article 

• In Section 2, we present the algebro-geometric background needed for our pur- 
poses. 

• We first review in Sections 3-4 the construction of an integrable system associated 
to a given compact Riemann surface. Its times are the moduli of a meromorphic 
1-form defined on the curve with fixed cycle integrals. The conserved quantities 
are the moduli of the curve. This construction dates back to Krichever [Kri77]. It 
encompasses the finite gap approach when the spectral curve is hyperelliptic, and 
produces in general the so-called algebro-geometric solutions of KP. We slightly 
reformulate this construction by defining a "spinor kernel" V'cK^i! ^2), globally 
defined on the spectral curve (Section 3.1). Baker- Akhiezer functions and the 
integrable kernel can be recovered from tp^i- We recall the construction of the 
Tau function %i and its properties. This integrable system is called " classical" , 
or " dispersionless" . 

• The goal of this article is to propose explicit formulas for another type of inte- 
grable system, called "dispersive" [Dub96]. The times are now (not necessarily 
independent) moduli of a spectral curve, which evolve at speed of order 
around a given algebraic spectral curve S. We define a function T which is an ad 
hoc solution to "loop equations" (Section 5), a spinor kernel ip{zi, Z2) via a Sato- 
like formula (Section 5.3), and correlators Wn{zi, . . . , Zn) encoding variations of 
T with respect to all possible moduli (Section 6). We present the conjecture that 
ip is self-replicating, and see its consequences: Hirota equations for T (Section 7), 
determinantal formulas for the correlators, Christoffel-Darboux formula for the 
integrable kernel, and a Lax pair formulation (Section 8). 

• Eventually, in Section 9, we recall how the results known about asymptotics of 
matrix integrals are compatible with special cases of our construction. 

2 Geometry of the spectral curve 

See [FK07, Fay70] for an introduction to algebraic geometry. 
2.1 Definition 

For us, a spectral curve S is the data of a complex curve C and two meromorphic 
functions X,Y on C: 

S = {C,X,Y) (2.1) 
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An essential example is given by the zero locus of a polynomial: 

8{x,y) = {) (2.2) 

If degy S = d, then for every value x, there are d values (with multiplicities) 
y = yi{x), . . . , yd{x) solutions of that equation. Those d functions can also be viewed to- 
gether, as d branches of one multivalued'^ function x H- y{x). The equation S{x, y) = 
can also be seen as the embedding of a Riemann surface C into CxCUjoo}. In other 
words, we prefer to represent the spectral curve in a parametric form: 

5 = {(x, y) e C X C I £{x, y) = 0} = {{X{z), F(^)) , z e C} (2.3) 

This means that there are two analytical functions (meromorphic) X and Y on C, 
which generate S. 

In this article, we only work with spectral curves for which C is compact. Then, 
there exists a polynomial S such that S{X,Y) = 0. We also assume C connected 
without losing generality. 

2.2 Some notations and properties 
2.2.1 Topology and holomorphic 1-forms 

The curve C is either simply connected, and then this is the Riemann sphere C = = 
C U {oo}, or it has genus > 0. Then, any maximal open contractible subset of C 
is called a fundamental domain. If it is of genus g > 0, there exist 20 independent 
non-contractible cycles (see Fig. 1), and we can choose them in such a way (but not 
unique) that: 

A n Bj = 6ij, A n Aj = 0, BiH Bj = (2.4) 

A basis satisfying these intersection relations is called " symplectic" . 

From the topological point of view, a genus > compact Riemann surface with 
a basis {Ai,Bi)i<i<g is a 40 closed polygon F, with edges 

[^1, B^, A^\ B^\ . . . , A, Sg, ^g"'] (2.5) 

glued by pairs. F is a fundamental domain of C. It is a classical result, that on a curve 
of genus 0, there exists independent holomorphic 1-forms dwj (holomorphic means in 
particular having no poles), and they can be normalized on the ^-cycles: 

j) duj = Sij (2.6) 

■*We could also consider the multivalued function y i— >■ x{y), and indeed, symplectic invariance 
guarantees that this would give the same Tau function. 
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Figure 1: A symplectic basis of 2q non-contractible cycles on a Riemann surface of 
genus g. 



Then, the x g matrix Tij: 

Tij = (p duj (2.7) 
is known to be symmetric Tij = Tj^i and its imaginary part is definite positive: 

r* = r, lmr>0 (2.8) 
r is called the Riemann matrix of periods of C. 
2.2.2 Riemann bilinear identity 

It is well-known that, for any meromorphic 1-form u defined on C: 



V Res u = 0. (2.9) 



P 

p pole of Lu 



This can be generalized to some forms which are not globally defined on C. Let Co'i,Co'2 
be two meromorphic 1-form defined on C, such that Res = for every pole of 1^2, 
and let 02 a primitive of UJ2 defined on a fundamental domain of C, i.e. d02 = ^2- For 
any symplectic basis of homology {A^ B) , we have the Riemann bilinear identity: 



^ Res 02u;i = ^ E (/ -1 ■ / C.2 - / . j 



UJ2 ■ (2.10) 



p pole of {(f>2i^l) 

2.2.3 Theta functions 

Given any symmetric matrix r such that Imr > 0, one can define the Riemann Theta 
function: 

0{u\t) = J2 e^^™-"e^" (2.11) 
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Since Imr > 0, it is a well-defined convergent series for all u in C^. Most often we will 
not write the r dependence of 9: 



e{\i\T) = e{u) (2.12) 

The Theta function has nice quasi-periodicity properties. If n, m G we have: 

^(u + n + rm) = e"*" ■"'■"W r-m) ^^^^ 13) 

It also satisfies the heat equation: 

d^. . 6 = — du-du (2.14) 
In this equation, Tjj and considered independent. 
2.2.4 Jacobian and Abel map 

Let us choose a generic basepoint o G C, which will in fact play no role. For any point 
z G C, we define: 

V2 G {l,...,g}, Ui{z) = [ du, (2.15) 



where the integration path is chosen such that it does not intersect any ^-cycle or 
i3-cycle. Then we define the vector: 

u{z) = [ui{z)]i<i<, GO (2.16) 

The application z t— u.{z) mod (Z^ + rZ^) is well-defined and analytical, it maps the 
spectral curve into the Jacobian J = C^/ (Z^ + rZ^). This defines the Abel map: 

z ^ u{z) mod (Z0 © rZ9) ^^'^^"^ 

The Jacobi inversion theorem states that every w G J can be represented as w = 
Sj=i for some points pi, . . . ,Pg E C. 

The Theta function can be used with r the Riemann matrix of periods of a Riemann 
surface C, and u the Abel map of a point on C. In this case, it enjoys other important 
properties. Its zero locus has the following description: there exists k G C^, so that 
6{w\t) = iff there exists g — 1 points zi, . . . , z^^i G C satisfying w = Yl%i ^i^j) + 
k is called a " Riemann vector of constants" , and it depends on the basepoint o used 
to define the Abel map u. 
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2.2.5 Prime form 



An odd characteristics c is a vector of the form: 

n + rm 



2 ' n,mGZ0, n*-mG2Z + l (2.18) 

The Theta function vanishes at odd characteristics: 6{c) = 0, and the following holo- 
morphic form: 

5 

1=1 

has only double zeroes on C, so that we can define its squareroot, and thus one can 
define the prime form: 

There exists choices of c such that E is not identically (we say c is "non singular"), 
and E is in fact independent of such c. It is a ( — 1/2, — l/2)-form on C, and it vanishes 
only at Zi = z^- In any local coordinate ^i^z) we have: 

i^(^i,^2) = 4=== +0((e(;^i)-e(^2))=') (2.21) 

zx~^Z2 ^d^(Zl) di{z2) 

Because of the Theta function, E{zi,Z2) is not a global (—1/2,-1/2) form on C. It 
transforms according to Eqn. 2.13. 

The Theta function associated to a Riemann surface satisfies a non-linear relation 
called Fay identity [Fay70]: for any Zi, Z2, z^, Z4 G C, any w G C^, 

^(w + c)^(ui2 + U34 + w + c) -^(^1' ^3)^(^2, ^4) 1 

E{zi, Z4)E{Z2,Z3) E{zi, Z2)E{zs,Z4) 

_ 9{w + U12 + c) 6{w + U34 + c) 6(w + ui4 + c) 6{w + U32 + c) 

^(2:1, 2:2) E{Z3,Z4) E{zi,Z4) E{Z3,Z2) 

where uji = u{zj) — u{zi). 
2.2.6 Bergman kernel 

The Bergman kernel [BS53] is the "fundamental (l,l)-form of the second kind" [FayTO], 
defined as: 

B{zi, Z2) = d,,d,2 In {e{u{zi) - u{z2) + c)) (2.23) 

It is independent of the choice of a non-singular, odd characteristics c. It is a globally 
defined, symmetric (l,l)-form, having a double pole at Zi = Z2 with no residue, and no 
other pole. It is normalized so that: 

B{-,z) = 0, f B{-,z) = 2i7Tdui{z) (2.24) 



(2.22) 
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Near zi = Z2, it behaves, in any local coordinate C,{z), like: 



B{zi,Z2) = j-T— 77r^+0{l) (2.25) 



We also define the fundamental 1-form of the third kind: 



r B{-,z) (2.26) 

J 22 



where the integration contour is chosen so that it does not intersect any ^-cycle or 
S-cycle. It is a 1-form in the variable 2, and a function of the variable zi,Z2., and it 
satisfies: 

/ d^,,,,, = , / d^,,,,, = 2i'K{uj{z^) - Uj{z2)) (2.27) 

It has a simple pole at z = zi with residue +1, a simple pole at 2; = 2:2 with residue 
— 1, and no other pole. In other words, in any local coordinate ^(2): 

dS.,,.M - f ^f, , (2.28) 

z^z^ ^{z) - ^{Zi) 

dS.,,.M - '"^^['-^ , (2.29) 

Notice that in the variable z it is globally defined for 2; G C (it has no monodromy if z 
goes around a non-contractible cycle), whereas in the variable zi (resp. 22) it is defined 
only on the fundamental domain, it has monodromies when 21 (resp. 22) goes around 
a non-contractible cycle Bj: 

dS^^+Bj,z2{z) = dSzi,z2{z) + 2mduj{z), (2.30) 

resp. 

d5'^i,z2+Bj(^) = dS:,^^z^{z) - 2iTiduj{z). (2.31) 

2.3 Parametrization of meromorphic forms 

2.3.1 Sheets, ramification and branchpoints, local coordinate patches 

If degy £ = d then £{x, y) = has d solutions, i.e., for every value x, there are d points 
2^(x), . . . , z'^{x) on the curve C such that: 

X{z\x)) = X (2.32) 

X is then a meromorphic function on C of degree d. 2'(x) is sometimes called the 
preimage of x in the i-th sheet. 
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Definition 2.1 We call "ramification points of order k", the zeroes of order k > 1 of 
the meromorphic 1-form dX . If a E C is a ramification point, the corresponding value 
X{a) is called a branchpoint. All the other points z E C at which X{z) is analytical, 
are called " regular points" . 

Definition 2.2 We say that a branchpoint Xa is simple if X~^{{xa}) consists in d — 1 
points, one of them being a ramification point of order 1, and all the remaining ones 
being regular points. 

2.3.2 Definition of local ccordinates 

Near a ramification point a of order k, ^ai^) = {X — X{a)Y^^''^^'^ defines a local 
coordinate on the curve. Simple branchpoints play a special role in Sections 3.7, 5.1 
and 6.1. For a simple branchpoint we have 

U^) = ^X{z)-X{a). (2.33) 

Since X is a meromorphic function of degree d, it has d poles with multiplicities, i.e.: 

oo^^...,oo^^ ^rfoo, = ci (2.34) 

i 

Near ooj, a good local variable is: 

UA^)=Xiz)-'/''-^ (2.35) 

Besides, we will need to consider also poles of a meromorphic form u. If is a pole of 
u, but not a pole of X, neither a zero of dX, a good local variable is: 

U^)=X{z)-X{p) (2.36) 

In this case, the multiplicity of p is dp = —1. We shall now always use the local 
coordinates ^{z) defined above. Notice that they depend only on the function X{z) 
and not on Y{z). 

Definition 2.3 Given a meromorphic 1-form U!{z) which has no pole at ramification 
points, let us call: 

V = {poles of u} 
'Poo = {poles of X} 

V = VUVoo 

To any p G P, we have associated a coordinate patch C,p on C centered on p. 
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2.3.3 Poles and times, filling fractions 

Definition 2.4 For any p G V, we define the "times" near p as the coefficient of the 
negative part of the Laurent series expansion of u near p: 



UJ[Z) = 

j>0 



I.e. 



tpj = Res u(z) ^pizV 

We also write collectively: 

Notice that the times tp Q = ReSp u are not independent, because the sum of residues 
of uj must vanish: 

$^W = (2.37) 

pev 

Definition 2.5 We define the "filling fractions" (also called "conserved quantities") 
by: 

2.3.4 Form-cycle duality 

To each time tk, we associate a differential meromorphic form Uk{z), as well as a dual 
cycle ul, and a dual orthogonal cycle ul-^: 

tk O Ukiz) ^ ul ^ ul^ (2.38) 

In such a way that: 

du{z) 



= ujk{z), uj{z) = '^tkUJk{z) (2.39) 

^(^) k 

Ukiz)= I B{-,z), tk= I CO, u*nu*^ = Sij (2.40) 

The symbol \x{z) Cleans that we differentiate keeping the local coordinates ^ fixed (their 
definition depend only on the function X). More explicitly we have: 

• Filling fractions ej — > ujj = 1st kind differential = 2i7rduj: 
ojj{z) = 2i7iduj{z) = B{ 
UJ* = B,, 
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Residues tp^ — > = third kind differential: 

ujpfl{z) = dSp^o{z) = B{z, 



^1 11 



^Cp, tp,o= ReSp(FdX), 



P,0 ~ 2ivP^ 

where o is an arbitrary basepoint on C, and Cp is a small circle surrounding p 
with index 1. As we mentioned, the tp^Q are not independent variables, and only 
,o)p^po ^ fixed pq are independent. As a consequence, we see that only 
differences Upfl — Upifl and uj* q — uj*, ^ are independent of a choice of basepoint o. 

• Higher times tpj with j > 1 — > oopj = second kind differential: 

ujpj{z) = Bpj{z) = Res z'^p^p{z')-^ B{z' , z), 
, ,* n 

/ — J_«'+i n 

The form a; is a linear combination of those meromorphic forms, and almost by defini- 
tion we have: 

Theorem 2.1 

s 

uj{z) = y^JkUJkjz) = ^2i7ieidui{z) + y^Jp,o dSp^o{z) + ^ tpj Bpj{z) 

k i=l pfzp p&V,j>l 

We also define its dual cycle 



2.4 The prepotential 

The fact that J * J * B(z, z') is symmetric, implies that there exists a function -Fo(t) 
called the "prepotential" such that: 

" " ^ c., 1^ " = " /■ ! B (2.41) 



'-^^i Jul* lyviiyvn J ^* J ^* 

The problem (this is why we write quotation marks) is that those integrals are not 
well-defined for times associated to third kind differentials. Such a statement is correct 
after an appropriate regularization. 

When z is in the vicinity of a pole p, we define: 
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It is such that uj — dVp has at most a simple pole at p. Given an arbitrary base point 
o G C, the following integral is well-defined: 



Lij{z) — (Wp{z) — t 



Vp{o) - tpfl\n^p{o) 



(2.43) 



/ip depends on the base point o, but only by an additive constant independent of p. 
Since ^ptp,o = 0, the sum Ylip'^pfil^p is thus independent of o. In some sense, Hp is 
a regularized version of j\ u (which does not exists). Since for all the other cycles, 
J^* UJ is well-defined, we can now define the prepotential. 

Definition 2.6 The prepotential is: 



1 



Res VpU + y^i ^p>o /Wp + f 



_p&V p&V 

Theorem 2.2 The first derivatives of Fq are given by, for j > 1: 



OF. 



and: 



UJ = Res ^ ^ u, 
p t' 

dFo^ 



dFo dFo 



dtpfi dtpi^Q 



UJ 



jjjp jjjpi 



The proof of this theorem has appeared in many works, in particular [Dub92, Ber03, 
?, Ber06, EO07a]. The definition of Fq is thus equivalent to write 

dF. 



k 

which means that Fq is homogeneous of degree 2. Another classical result is: 
Theorem 2.3 The second derivatives of Fq are given by: 

d'FQ 



(2.44) 



except for the following cases: 



d 



d 



dtkdti 



d 



B 



9tk \dtpQ dtpi 



d 



d 



d 


d 


dtpfi 


dtpifi 


' d 


d 


^dtpfi 


dtpi^Q 



dtp^o dtpi^Q 
d d 



i^p,0 Otp"fl 

d d 



dtpfl dtp/. 



^0 

Fq 
Fq 
Fq 



B 



k p,0 



B 



In {Eip^p'fdUp) d^p'ip')) 



In 
In 



E{p,p')Eip,f) dUp) 

E{p',p") 
E{p,p')E{p',p) 
E{p,p)E{p',p') 
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The second derivatives of Fq do not depend on the form u, and thus do not depend on 
the times. Thus we have: 

P.45, 

3 Reconstruction formula 

In this section, we review the reconstruction of a Lax matrix and an isospectral inte- 
grable system from the spectral curve. This so-called algebro-geometric reconstruction 
was first discovered by Krichever [Kri77], and we merely rephrase it. The only differ- 
ence is that, instead of Baker- Akhiezer functions, we prefer to use a " Baker- Akhiezer 
spinor kernel", which turns out to be a more intrinsic object. 

3.1 Baker- Akhiezer spinor kernel 

Define the 1-form: 

t) = ^p.o dS'p,o(z) + 5^ XI 
p€V pev i>i 

= uj{z) — 2i7r 'Y^eidui{z) (3.1) 

i=l 

which depends linearly on the times (and not on the filling fractions). By construction 
X is normalized on ^-cycles: 

<f X = (3.2) 
-'A 

Then we define the vector ({t) = [Cj(^)]i<i<g with coordinates: 

^'CT = ^ni^ = ^n(i^-t i -) (3-3) 
which we write for short as: 

at) = ^ / (3.4) 

It can be decomposed as: 



Uk (3.5) 



pev i>o k={p,j) 



The vector ({t) is a linear function of the times and is independent of the filling 
fractions e^. In other words, it follows a linear motion with constant velocity in 
the Jacobian, as a function of any of the times. A well-known property of integrable 
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systems is that, in appropriate variables, the motion (with any of the time tk) is uniform 
and hnear. The algebraic reconstruction takes the linear evolution in the Jacobian of C 
as starting point, and produces more intricated quantities whose evolution is described 
by a Lax system. 

Definition 3.1 We now define the Baker- Akhiezer spinor kernel as the (1/2,1/2) 
form: 

^ ' E{z„z,) e{at) + c) 

where c is a non-singular, odd characteristics. 

We write a subscript ci to distinguish the algebraic construction (associated to a "clas- 
sical" or dispersionless integrable system) from the one proposed in the second part 
of the article. This kernel is also known as the Szego kernel [Sze21], twisted by the 
exponential term e-^^a 

Theorem 3.1 il)c\{zi,Z2] i) is a globally defined spinor in [zi,Z2) E C xC, i.e. it is 
the squareroot of a symmetric (1, l)-form. 

• It has a simple pole at Zi = Z2 : in any local coordinate ^ (z) 



i>ci{zi, Z2 ; t) - - ■ ^ 



zi^Z2 E{zi,Z2) ^{Zl)-^{Z2) 



• It has essential singularities when Z\ (resp. Z2) approaches a pole of u. 

Proof. The behavior at zi — )■ Z2 is obvious, and the essential singularities at the poles 
of u come from the exponential term. Next, we need to prove that ipci{zi,Z2; i) is 
unchanged when zi (resp. Z2) goes around a non-trivial cycle. When zi (resp. Z2) goes 
around an ^-cycle, the vector u{zi) (resp. u{z2)) is translated by an integer vector, 
9 is thus unchanged, and thanks to Eqn. 3.2, ipci is unchanged when zi (resp. Z2) 
goes around an ^-cycle. When zi (resp. Z2) goes around a i3-cycle, the vector u(2;i) 
(resp. u.{z2)) is translated by a lattice vector of the form r ■ n with n G Z^, and 6 gets 
multiplied by a phase according to Eqn. 2.13. Remember that the prime form E{zi, Z2) 
is also a 9 function, and also gets a phase given by Eqn. 2.13. ipd is thus changed by: 

Mzi + nS, Z2;t)^ Mzi. Z2 ; t) CC*") e^'hx (3.5) 

and because of Eqn. 3.3, i.e. C = ^ j^X^ see that tpd is unchanged when zi (resp. 
Z2) goes around a S-cycle. □ 
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3.2 Duality equation 

Then we construct the following spinor matrix of size d x d: 

^ci(a:i,X2; t) = [M^'{xi),z'{x2);P)]ij=, (3.7) 

where we recall that z'^{x) are the d preimages of x on the curve C, i.e. X{z'^{x)) = x, 
and d = degX. These preimages are distinct and this matrix is well-defined when xi 
(or X2) is not at a branchpoint. 

Theorem 3.2 We have the "duality" equation: 

Wci(xi, X2 ; t) Wci(a;2, ; t) = W<.i(xi, X3 ; t) 

(Xi - X2)(X2 - X3) 

Proof. ^ 

dX{z) ^^^(^'(^i)'^' t)?/;ci(2;, 2;^(x3) ; t) (3.8) 

is a meromorphic function of z. Indeed, the product of two (l/2)-forms is a 1-form, 
and when we divide by dX, we get a function. The essential singularities coming from 
the exponentials cancel in the product, so this function can only have poles, i.e. it is 
meromorphic. The only possible poles are at 2; = z^{xi) or z = z^{x3) or at the zeroes 
of dX{z). Then, summing over all sheets, we see that 

E ipdjz'ixi), z^{x2) ; t) ipciiz''{x2), z^jxs) ; i) 
^ dX(2fc(x2)) ^ ' 

is a symmetric sum of a meromorphic function over all sheets of X2, therefore it is a 
meromorphic function of X2 G C, i.e. a rational function of the complex variable X2- It 
remains to find its poles. l/dX(2^(x2))) behaves like 0(x2 — X(aj))^^/^ at ramification 
points, and since a rational function of X2 cannot have a singularity of power —1/2, this 
means that this rational function has no pole at branchpoints. Its only poles can then 
be at X2 = Xi or X2 = x^, and they are simple poles. The residues of the corresponding 
poles are easily computed and give the theorem. □ 

Theorem 3.3 We have a refined version of the duality equation: 

iJciizi, z ; i) ipci{z, Z2]t) = -il)c\{zi, Z2 ; t) {^S^^^^^{z) - 2i-n ^ aj{zi, Z2 ; t) duj{z) 
where 

e.^^ (u(2;i) - ^x{z2) + C(t) + c) Qu, (C(i") + c) 

aj{zx,Z2 ; t) 



^(U(zi) - V.{Z2) + C(t) + c) ^(C(t) + c) 
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Notice that Theorem 3.2 is a corollary of Theorem 3.3. Indeed the duality equation 
(Theorem 3.2) can be obtained by summing the equation above on all sheets z = z''{x), 
because '^i.dui{z^{x)) = and 

{X{z,)-X{z2)) dX{z 



Proof. Notice that ijjci{zi^z; t)ijjc\{z, Z2', t) is a meromorphic 1-form in z, since it 
has no essential singularity. It has simple poles a.t z = zi and z = Z2, with residues 
^ipciizi, Z2; t), and it has no other pole. This means that ipciizi^z; tjipdiz, Z2:, i) + 
■^01(^1,^2; t) dSz^^z2iz) is a holomorphic 1-form, with no poles, therefore it must be a 
linear combination of the dui{zys, which we choose to write: 

V^cK^i, z ; t) il^diz, Z2]t) = - ipc\{zi, Z2 ; t) {^Sz^^z^{z) - 2i7i ^ aj{zi, Z2 ; i) duj{z) 

(3.11) 

The left hand side is a well-defined spinor of Zi and Z2, whereas in the right hand side, 
dSzj^^z2{z) = f^^B{z,-) gets some shifts when zi or Z2 go around non-trivial cycles. 
This implies the following relations for the coefficients aj{zi, Z2): 

aj{zi+ Ak,Z2; = aj{zi,Z2; i) , aj{zi, Z2 + Ak ; i) = ai{zi, Z2 ] i) (3.12) 

aj{zi + Bk,Z2] i) = aj (21 , 22 ; i) - 21716 j^k , aj{zi,Z2 + Bk] i) = aj {zi,Z2] i) + 2i7r6j^k 

(3.13) 

Moreover, we must have aj{zi,Zi; t) = 0, and aj{zi,Z2:, i) may have poles when 
'4>ci{zi,Z2; i) = 0. Apart from those poles, aj{zi,Z2:, i) has no other singularities. 
The following quantity has all the required properties: 

9^^{u{z^)-u{z2) + C + c) 9^^{C + c) 

^(u(zi)-u(z2) + C + c) eic + c) ^ - ' 

So, the difference of aj and that quantity should be a meromorphic function of zi and 
Z2 without poles, i.e. a constant, and its value is zero by looking at zi = Z2- □ 

3.3 Link with Baker- Akhiezer functions 

Given g points pi, . . . ,pg G C in generic position, the Baker- Akhiezer function ipBA^z:) 
is usually defined as the unique (up to a multiplicative constant) function on C, with 
poles at the points pi, . . . ,Pq, and essential singularities at the poles g of a meromorphic 
1-form cr on C, such that In '^^^{z) = cr+0(l) when z ^ q. It can be readily written 
as: 

Wba = exp / a) — — — — 3.15 
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We can obtain such a function by specializing the spinor kernel associated to the 1-form 
cr: 

^BAl^) = , (3.16) 

^<lX{z)<lX{p,) 

Indeed, by Jacobi inversion theorem, we can always represent u(pi)— c— k as X]^=i ^i^j) 
for some points gi, . . . , gg G C. Moreover, according to the description of the zero set 
of the Theta function, gi, . . . , qj are the q zeroes of 6{yi{z) — \i{zi) + c), in particular 
Pi e {gi,...,g0}. 

3.4 Eigenvectors and Baker- Akhiezer functions 

The usual formulation of integrable systems is obtained by specializing one of the points 
to X = oo. In some sense, we would like to consider: 

^cIk(^)" ='>ci(^,oo,) (3.17) 

The problem is, that the expression in the right hand side is divergent, and thus we 
again need regularization. 

The definitions in this paragraph also apply to the spinor kernel constructed in 
Section 5.3, so we drop here the ci index. We define: 

V^,o(z)= hm ^%^y-^^'^HiooX^2)f--^ (3.18) 

and if c/oo, > 1, we define for < j < (doo, — 1): 



(3.19) 



There are d pairs / = («, j) such that < j < d^oi — 1, and therefore the vector: 

V^(z) = [i,i{z)] (3.20) 
is a d-dimensional vector, and the matrix: 

^{x-^ = [iji{z\x))]i,i<k<d (3.21) 

is a d X d square matrix. 

Similarly, we define the dual Baker- Akhiezer functions: 

4>^,{z)= hm 4=4^«"''°°^^"^(^-(^i))-*-'° (3.22) 
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and if doo, > 1, we define for each < j < (doo, — 1): 



(3.23) 



Zl=COi 



There are d pairs I = such that < j < di — 1, and therefore the vector: 

0(z) = IM^)] (3.24) 
is a (i- dimensional vector, and the matrix: 

^x;^) = [^i{z'{x))]i,i<k<d (3.25) 

is a d X d square matrix. 

3.5 ChristofFel-Darboux relations 

Theorem 3.4 The matrix 



dx 



is invertible, and independent of x. The matrix A^i is called the Christoffel-Darboux 
matrix. This can also be written: 

**i(x) Ad$ci(x) = id dx 

Proof. This is a mere apphcation of Theorem 3.2, up to a conjugation. Indeed: 

{A~i^)(i,k),{i',k') 



■*oo./,0 



Zl=00-,,Z2=OOi 



de^;n^i) de-n^2) L \X{z) - X{z,) X{z) - X{z,) 

MZUZ2) e^°°»(^^)-^°°.'(-^) eoo.(^2)*°°-°eoo,(^l)"'°°'"° - 

^d^oo,,(2;i) d^ooi(2;2) Jzi=oOi,,^2=oOi 

If 2 7^ i', the following quantity 

Mzi,z2) e^°°'^^^)-^-.'(^^) eoo.(^2)*°°-°eoo,(^i)"'°°'"° 

Vd^oo,,(2;i)d^oo,(2;2) 



(3.26) 



(3.27) 



19 



has a well-defined limit when zi — )■ ooi' and Z2 — )■ ooj, and the term 



1 



1 



X{z)-X{z^) 



X{z)~X{z2) 



behaves like: 



(3.28) 



X{z) — X{zi) 2i-s>oo., 

so we are computing the [k' — l)-th derivative of 0{^oo^X^iY'^^')i where k' < doo^,, and 
therefore we get 0, i.e.: 



cl\ii,k),{i',k') 



iii^i' 



li i = i', we first take the limit zi — ooj, and again the term with 



(3.29) 
vanishes. 



X{z)-X{zi) 

Then, remember that ipci{zi,Z2) has a simple pole at zi = Z2, and thus the derivative 
with respect to zi, can generate a pole of degree k' at Z2 = oOj. Therefore, we are 
computing the {k — l)-th derivative of 0{^ocX^2Y°°'~''')- We get zero k + k' < doo^, 
and therefore: 

(^"')cl|{^,fc),(^,fc') =0 ifk + k'<d^^ (3.30) 

If i = i' and k + k' = doa + 1, the only non- vanishing contribution is: 



— 7— TT — T lim 

d^i). (22) ^1^°°- 



(U(^i)-eoo.(^2))'^' 

Mzi,Z2) E{zi,Z2) e^-«(^^)-^-«(^^) UX^2y-'-" UX^i)-'--' 
lim ilj{zi,Z2) E{zi,Z2) ^.{^2)'^-' UXzi)' 

= {-lf'{k'-l)\{k-l)\ 

^ 

The matrix A^-^ has thus typically the shape: 

/ 

* * * 



Z2 = OOi 



(3.31) 



^cl 



v 



(3.32) 



It is made of (Inverted) triangular blocks. Since the diagonal of each triangle is non- 
zero, this proves that the matrix is invertible. 
Then, ii i = i' and k + k' > d^o + 1, we write that: 



X{z)-X{z,) Xiz,) 



Oil/Xiz^f 



(3.33) 
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and we see that the leading term X(I7) independent of X{z), and the part which 
depends on X{z) is 0{l/X{ziY) = O(,^ooi(^i)^'^°°0- ^ ^'^^ vanishing contribution to 
the part which depends on X{z) could occur only ii k + k' > 2dao,, which can never 
happen since we assumed k, k' < doa- This proves that is independent of X{z). □ 

Corollary 3.1 The matrices \E'ci(a^;0 ^'^^ ^ci{x',t) are invertible. 

As a consequence, ipci{zi,Z2) can be identified with an integrable kernel. This means 
that we have: 

Theorem 3.5 Christoffel-Darboux relation: 

E / J ^cl|/(2:i) Acl\IJ 0cl| 7(2:2) 



^01(^1,2:2) 



X{Z,)-X{Z2) 

Proof. This is an application of Theorem 3.4 and Theorem 3.2. Indeed, we have from 
the duality equation, for any x, and for any invertible matrix C: 

{X{zi)-X{z2))^ci{X{zi),X{z2)) 
= ^,iXiz,U)^^ix,Xiz2)) 

= <iJ,yiXiz,),x)C{x) C-\x)^,ix,X{z2)) (^(^1) -X) 

dx 

The very definition of the ipci\i^s, means exactly that there exists a matrix Cc\{x) such 
that: 

^*i(X(zi)) = lim Cci(x) (3.35) 

and similarly, there exists a matrix Cd such that: 

^ci{X{z2)) = lim C,,{x) *ci(x,X(^2)) (3.36) 
When zi = Z2, we have \l/ci(a;, X(2i))\l/ci(X(2i), x) = — ^!l_'^(ll)yi id, which implies: 

- hm /"^^^Z'^^,, Cd{x) Cd{x) = $ci(X(zi)) ^*,(X(zi)) = dX(.i) (3.37) 

x^co \^X — A \Zi)j^ 

and therefore: 

C:,\x)=A,,C,,{x)dx (3.38) 

This means that: 

{X{Z^) - X{Z2))^,X{X{Z^),X{Z2)) 
= lim V^el(X(^i),x)Cel(x)AelC7el(x)vl>,i(x,X(^2)) (^(^^ " ^)(^(^2) - x) 

= ^*i(X(^i))Ai$ci(X(^2)). (3.39) 

□ 
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3.6 Lax matrix 

As corollary of Theorem 3.2: 

Corollary 3.2 If Xi and X2 are not branchpoints, and if X2 7^ Xi, then the matrix 
\I^ci(3;i, X2) is invertible, and: 

^ T / -\ dx2 

Wci(a;i, X2 ; t)Wci(a;2, xi;t) = - id 

(xi - X2y 

Proof. Take 0:3 = xi in the duality relation. □ 
Corollary 3.3 Reconstruction formula. Let 

L{x) = diag(n(x), . . . , = diag(F(^^(x)), Y{z\x)), . . . , F (/(x))) 

For every xi, the matrix: 

"I 2/")^ 

L^i (x; i) = ^ci(a;i, x ; t) L{x) '^~^{xi,x; t) = — - — ^'ci(a;i, x; i) L{x) ^'ci(x, xi, i) 

dxi ax 

is a rational function of x. Its characteristic polynomial is independent of the times t, 
and its zero locus defines the spectral curve: 

det [y - L^^ (x; t)) = det{y - L{x)) = £ (x, y) 

Changing Xi just amounts to a conjugation: 

Lx'^{x]t) = ^'ciK,a;i;t) L^^{x]t) ^^^(x'l, Xi; t) 

Theorem 3.6 The matrix Lx^{x;t} obeys the Lax equation: 

d 

— L^, (x; t) = [Mt^.^^ (x; t) , L^, (x; t)] 

with respect to any of the times tk, and where the matrix Mti^-xi is: 

d 

Mt,,x^{x;t) = — ^'ci(a;i,x;t) ^'^-J^(xi, x; t) 
dtk 

Therefore, from any plane curve C : S{x,y) = 0, we have constructed a Lax system. 
Since the spectral curve here coincide with C and is independent of the times, the 
evolution is said " isospectral" . This reconstruction of the Lax matrix from the spec- 
tral curve, is well-known in the theory of integrable systems [BBT02]. Usually, it is 
presented with special values of xi, for instance xi is most often chosen as a pole of X. 
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3.7 Differential systems 

Baker- Akhiezer functions satisfy simultaneously several first order differential systems 
with respect to the spectral parameter and the times. 

Theorem 3.7 The matrix \l/ci(a;i, x; t) is solution of a linear ODE with respect to the 
spectral parameter x: 



dx X — Xi ' / ^/dxdx 

where M^-^^i is a rational function of x having poles only on V , with the same order as 
those of uj{z) /dX{z) . 

Proof. Write x = X{z), we have: 

dX{z) 



dXiz) [M,,,,iXiz);^]^^^^-6,, 



{X{z) - xi) 



{x{z) - x,r V d I --7:;::^: ^ i "7:j::7" (3-40) 



k 



^/dX{z) dxi I ^/dX{z) dxi 



In the right hand side, the essential singularities cancel, and only meromorphic singu- 
larities remain. Since we perform the sum over all sheets, the result is necessarily a 
rational function of X{z). Poles could occur at singularities of ip, or also at zeroes of 
dX{z), or at X{z) = Xi. 

ipci has a simple pole a.t x = xi, i.e. at z^ = z^{xi) for some k. Taking the 
derivative yields a double pole, and multiplying by {X{z) — xi)^ cancels the double 
pole. If z = j, there is also a simple pole coming from ipci{z^, z^{xi)), and if z 7^ j there 
is no pole. It is easy to check that the residue of the pole at a; = xi is —Sij. This 
implies that M^-xiix) has no pole at x = Xi. From its definition, Mx^x-^{X{z)) must 
behave as o{{X(z) — X(a))~^/^) at a ramification point a, and since it is a rational 
fraction of X{z), this must actually be 0(1), meaning that M^.^-^lx) has no poles when 
X approaches a branchpoint. 

Near z ^ p & V, the only singularity comes from the exponential term and we have 

+0(1) (3.41) 

This shows that dx Mr^-rj.^[x; t) has poles in V, of order at most that already present in 
u/dX. □ 

Theorem 3.8 We have with respect to the times tpj (j > 1): 
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{{} {dt^^ - Mt^ .^,{x ; t}) ^ei(a;i, x;t) = 0. 

iii) Mtpj-xi is a rational function of x, with possible poles only at x = X{p). 
IfX{p) oo, the pole is of order j . IfX{p) = oo, the pole is of degree 1+ [{j — '^)/dp\ . 

One could write a similar theorem for dt^Q — dt^, ^. 
Proof. We have: 

= {M^'{xi),z''{x))) M^'{x),z^{x,)) (3.42) 

k 

and again the essential singularities cancel, and since we sum over all sheets, the result 
is a rational fraction of x. Since the pole of ipci at x = Xi has a constant residue, 
the derivative with respect to tpj has no pole at x = Xi. Since we divide by dX{z), 
poles at branchpoints could occur, but dX{z) vanishes as ^yx(z) — a, and a rational 
fraction of x cannot have half-integer singularities, so there is no pole at branchpoints. 
Therefore, poles may only occur through dt^.e^'^^^, and we find: 

[M,^^.x,{x-t)\^ = -^^^^ E (£(7"""") M^\x),z^{x^)) 

+0(1) (3.43) 

which shows that Mt^ .-^^lix'-it) has poles only at x = X(p), and of order at most j. If 
X{p) = oo, then (x — Xi)^/dx has a pole of order dp — 1, and thus Mt^ j-xi{x:,t} has 
a pole in z of order at most dp + j — 1, and thus in x it has a pole of order at most 
{dp + j ~ ^)/dp- Since Mt^j-xiix; i) is a rational fraction of x, the order of its pole must 
be integer and is thus at most 1 + [(j — l)/dp\. □ 
Again, in the usual presentation of integrable systems, the ODE's are presented 
after sending xi to oo, i.e. in terms of \l/(x; t) and $(x; t). 

3.8 Isomonodromic problem 

We see that the matrix '^cii^i, x;i) is a simultaneous solution of several compatible 
linear differential systems with rational coefficients. This leads to some isomonodromic 
problem. 

We have 

d 

— "^ciixux;^ = M^;^.,(x)^'ci(xi,x;t) (3.44) 
ox 

A general property of linear differential systems, is that the solution '^^lixi, x;i) is 
analytical wherever M^-xi (x) is analytical, and may have essential singularities at sin- 
gularities of Mx^xiix). The solution '^ciixi, x;i) may have monodromies around the 
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singularities. Let p be one of the poles, then the monodromy after going around p is: 

^ci(a;i, t) = '^c\{xi, x_; t) Sp,^, (3.45) 

where x_ and x+ mean the same value of x before and after going around the singularity. 
Since Mx-xi {x] t) is a rational function of x, it has no monodromy and we have 

M,.,,(x+;r) = M,.,,(x_;t) (3.46) 

and it is easy to see that this implies that 

^^P,^.=0 (3.47) 

In other words, the monodromy is independent of x, this is why the system is called 
"isomonodromic" . Similarly, consider any time tpj-, and write: 

d 

^ci(a;i,a;;t) = ..^,(x) ^ci(a;i, x; t) (3.48) 



Since Mt^..^xx{x) is a rational function of x, it has no monodromy, and this implies that 

7^ = (3.49) 

In other words, the monodromy is independent of the times tp^j. Here actually, the 
monodromy group has finite order since \l/(xi,x;t) lives in a bundle over the fixed, 
compact algebraic curve C). This is a very special case, which is implied by the fact 
that the evolution is isospectral [BBT02]. 

4 Dispersionless Tau function 

Definition 4.1 We shift the prepotential and define: 

i=l * j,i'=i 
= Fo-C-e-^e-T-e 



1 d'Fo 1 

Definition 4.2 We define: 



2 ^ dtkti 2 



where u = t^Uk is the 1-form depending linearly on the times and the filling frac- 
tions, and Fq is the shifted prepotential associated to u. We shall also use as notation 

TciiC,u) = Tcii^ 
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depends as well on the data of a non-singular odd characteristics for the Theta 
function. It is the Tau function associated to the solution of the problem of isospectral 
evolution described in Section 3.6. We call here Tau function, any function satisfy- 
ing Hirota bilinear equation, and linked to Baker- Akhiezer functions by Sato relation 
[JMUSla, JMUSlb, JMUSlc]. Sections 4.1 and 4.3 prove that %i fits in this definition. 
In this article, we also call Td the " semiclassical" Tau function. 

4.1 Sato relation 

Sato relation [SS83] means that the Baker-Akhiezer kernel ipci{zi,Z2) is obtained from 
the Tau function by a Schlesinger transformation. This can be formulated as: 

Theorem 4.1 

Proof. We add two simple poles at zi and Z2 by considering: 

ux{z) = u{z) + \dS^^^Z2{z) 
According to Theorem 2.2 we have 



dX 



A=0 



U) 



Z2 



and according to Theorem 2.3 we have 

= -ln {E{zr,Z2fdX{zr)dX{z2)) 



A=0 



Since Fq is homogeneous of degree 2 we have 

Fo{ux) = Fo{u) + X 
= Foiu 

It is also easy to see that 

A 



+ A 


dFo 




A2 


d'Fo 


dX 


+ 

A=0 


2 


dX^ 


+ A 




In 


{e{zi 




J Z2 









A=0 



2in 



C ^ C + 77- f dS,,,,, =C + A(u(zi) - u{z2)) 



B-tA 



Taking A = 1 implies the theorem. 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 

□ 
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4.2 Expansion near poles 

In order to recover the usual presentation of Sato relation, let p & V, and Up C an 
open neighborhood of p on which the local coordinate C,p is well-defined. The Schlesinger 
transformation: 

u u + dSz^^z2 (4-6) 
can also be written by Taylor expanding dS^^^^j Zi, Z2 G Up-. 

= 2^ — — ujpj - Upj (4.7) 
j>i 

Thus, the times associated to p after Schlesinger transformation are: 

Vi > 1 + (^ - &^) (4.8) 
In the literature, a notation [z]p is introduced for this special collection of time: 

N,= (^l (4.9) 



\ J / j>l 

And since dSz^^z-zi^) has vanishing ^-cycle integrals, the filling fractions are not 
changed. For example. Theorem 4.1, together with Eqn. 3.18 gives: 

/ / N '7^1 + [^l]p) /A 1n^ 

'^ci\{p,o){zi) = —77 (4.10) 

lc\{t) 

4.3 Hirota bilinear equation 

For any point 2; G C, we define the insertion operator 6^, which acts on functions of a 
meromorphic 1-form on C, by considering the one parameter deformation consisting in 
adding a Bergman kernel B{z, ■)/dX{z): 

../(.)^dX(.)A;(.+ A|^)[^^ (4.11) 

Locally, if z is near a pole p, we have: 

6z is obviously a derivation. 

Theorem 4.2 For any two 1-forms u and Cj defined on the same Riemann surface C, 
we have: 

Res Vci(^, z' ; u) ^l)c\{z\ z; Cj) = -6^ In '^''^^^^ 
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Proof, ipdiz, z' ; u)ipci{z', z ; Cj) has a double pole at z = z', and thus evaluating the 
residue computes a derivative. □ 

As a corollary, if we choose u = u + dSz-i^,z2J the expression on which we take the 
residue is a meromorphic form in z (it has no essential singularity). Its only other pole 
is located a.t z = Z2, it is simple, and by moving the integration contour, we find: 

Corollary 4.1 For any data {C,u), the Baker- Akhiezer kernel is self-replicating: 

4 ^cl(^l, Z2) = -i^c\{zi, Z) i)c\{z, Z2) 

This self- replication property is the analog of the Ricatti equation in [BBT02]. Notice 
that it can also be obtained by a straightforward computation of 5z\ii-ipc\{zii Z2) from 
Definition 3.1) and comparison with the refined duality equation (Eqn. 3.3). 

Associated to any derivation dt, we may define a Hirota operator Dt [Hir71] acting 
on two functions f(t),g{t), such that: 



(A/ ■ gm ^ dJit + u)git -u) = g\t)dt{f/g) (4.13) 

■u=0 

This allows us to reformulate the self-replication property as: 

Theorem 4.3 In terms of%\, for any polep, the self-replication property is equivalent 
to the Hirota bilinear difference equation: 

Dzrcl{t+ [Z^]p - [Z2]p) ■ Tclii) = -%l{t+ [Zi]p - [Z\p)rcx{t+ [Z\p - [Z2]p) 

Actually, this property for Tc\ is equivalent to the Fay identity (see § 2.2.3) satisfied by 
the Theta function 0{-\t). The Hirota equation can also be written in a more symmetric 
way by setting: 

r, r , Nip 



(4.14) 



Namely: 



2 2 y V 2 2 
= -7;, [t- [z\, + Mi + ((-+ - M _ M 

The procedure to translate it into a set of differential equations (with respect to the 
times) is well-known, it is merely obtained by Taylor expansion in ^p(^), ^p{zi) and 
ip{.Z2) (see e.g. [JM83]). This gives an infinite set of partial differential equations 
involving derivatives of the Tau function with respect to times tpj, which are equations 
of a KP hierarchy. These equations are equivalent to Hirota bilinear equation. 
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5 Dispersive deformation 

So far, the spectral curve (C, X, Y) was fixed once for all, and was equipped with a 
1-form u = tkOJk depending linearly on times. Now, we shall let the spectral curve 
itself change around (C, X, Y), and in particular we may vary the complex structure of 
the curve C. YdX shall play the role of the 1-form. 

5.1 More geometry: symplectic invariants 

For any spectral curve S = (C, X,Y), a sequence of symplectic invariants Fg{S), and of 
symplectic covariant forms uJn\S) was defined in [EO07a]. Let us recall their definition 
and main properties. 

5.1.1 Topological recursion 

Let flj be the ramification points of S, i.e. the zeroes of dX. We assume that the 
spectral curve S is regular, i.e. are simple zeroes of dX and dY{ai) ^ 0. Then, Y 
behaves like a squareroot near ramification points: 



Y{z) = Yiai) + r (a,) ^/Xiz) - X(a,) + 0(X(^) - X{ai)) (5.1) 

and there is a unique other point z ^ z such that X{z) = X{z), at least for 2; in a 
neighborhood [/q- of aj. Then we define the recursion kernel [EO07a]: 

^^''^ - ~2iYiz)-Y{z))dX{z) ^^-^^ 

In the variable zq, it is a meromorphic 1-form defined globally for zq & C, and in the 
variable z, it is the inverse of a 1-form, only defined in [/q. . 

Definition 5.1 The "symplectic covariant forms" ujn\zi, . . . , Zn) are meromorphic 
(1, . . . , l)-forms defined by the following recursion: 

n times 

Jf\z) = Y{z)dX{z), a;f (^1,^2) = B{z,,Z2) 

and for J = {zi, 

UJnll {zo, J) = V Res K{zo, z) (uj^^^2 {z, z, J) 

i 

'^stable" 
0<h<g,ICJ 

where "stable" means that terms containing aoj^^ factor should be excluded of the sum. 
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An important property proved in [EO07a], is that: 

Theorem 5.1 For 2 — 2g — n < 0, uin\zi, . . . , Zn) is symmetric in its n variables, and 
it is a meromorphic form in each variable, having poles only at ramification points, 
with vanishing residues. The poles are of order at most 6g + 2n — 4. 

For instance, for {g,n) = (0,3), this definition amounts to: 
_ P B{zo,z)B{zi,z)B{z2,z) 

~ ^.^l dX(^)dF(^) ^^-^^ 

Definition 5.2 The "symplectic invariants" Fg are numbers associated to S, as fol- 
lows. 

• For g = 0, we define Fq{S) as the prepotential, see Definition 2.6. 

• For g = 1, Fi is defined in terms of the Bergman Tau function Tb,x introduced 
by Kokotov and Korotkin [KK03]: 



where Y\ai) = liuiz^ai — Y{ai))/ ^yX{z) — X{ai). Fi is related to the 

logarithm of the determinant of the Laplacian on C with metrics |?/dxp. 



• For g > 2: 

^ i 

where ^{z) is any primitive ofYdX, i.e. d$ = YdX . 
As a notation we write cJq^^ = Fg. 

The name "symplectic invariants" came from the important property, proved in 
[EO07b], that Fg and the cohomology class of Un^ are invariant under the follow- 
ing transformations of spectral curves, each of which let the symbolic symplectic form 
dX A dY invariant: 

(X, Y) ^ X), (X, Y) ^ (X, y + R(X)), (X, Y) ^ (AX, Y/X) (5.4) 
where R is any rational function and A G C. 
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5.1.2 Special geometry 



Let us consider an infinitesimal deformation of the spectral curve S. It implies an 
infinitesimal deformation of the complex structure of the Riemann surface C, and a 
deformation of the form YdX: 



YdX — > YdX + XQ + O(A^) 
where f2 is a meromorphic form. We can decompose'' as in Theorem 2.1 

r2 = tk cok + 2«7r ej duj 
k j=i 

There is a dual cycle associated to ^2: 

k j=i 



such that: 



In* 



(5.5) 



(5.6) 



(5.7) 



(5.^ 



Theorem 5.2 

Vra > 



d_ 



A=0 



and the derivative is evaluated at X{zi) fixed. 



This theorem was proved in [EO07a]. This set of relations is usually called "special 
geometry. Let us give some examples: 



1st kind deformations/ addition of holomorphic forms: 

duji^\zi, ...,Zn) 



dei 



OJ^^lli', Zi, . . . , Zn) 



(5.9) 



second kind deformations/addition of a double pole: 

If p is a pole of YdX, this does not apply to {g, n) = (0, 0). 



(5.10) 



^LOk here belongs to the basis of meromorphic forms defined in Section 2.3.4. It has nothing to do 



with the invariants 



,(9) 
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third kind deformations/ addition of simple poles: 



duji^\zi, ...,Zn) _ duJn\zi, . . . , Zn) _ 

dtp^o dtpifi Jp, 



V 



^nlA-iZl,...,Zn) (5.11) 



{g,n) = (0,0). We recover the constitutive relations for the prepotential (Sec- 
tion 2.4): 

dFo 



YdX (5.12) 



k 



{g, n) = (0, 1). We recover the definition of the form-cycle duality pairing by the 
Bergman kernel: 

r 

-iz)=Uk{z) = f Bi;z) (5.13) 



dYdX 



5.2 Candidate for a dispersive Tau function 

Let us consider: 

[S, n] = exp iV2-2, p^^^^ ^ ^ ^) j 

where Fg{S, n') are the symplectic invariants associated to the curve S and the vector 
of filling fractions n' = {n[, . . . ,n'^). lnZ[iS,n] is a formal Laurent series in N (in 
particular we emphasize that it contains no oscillatory terms, by opposition with the 
definition to come). Then, we would like to define the sum over all integer filling 
fractions: 

7[,,.](5)= 5^ e2-"-Z[,][5,n] (5.15) 

The intuition behind this definition is Whitham averaging. The filling fractions are 
considered as fast evolving moduli in the system, so that we do not see a fixed spectral 
curve, but rather a statistical mixture of spectral curves with different filling fractions. 
When > 0, In T[fj.^u]{<S) is no more a Laurent series in A^. To give a precise meaning 
to this expression, one has to collect the terms order by order in 1/A^. 

Initially, such an object has been introduced with heuristic arguments to compute 
the large asymptotic of hermitian matrix integrals, and it leads to theta functions 
to the leading order [BDEOO], and a whole series of corrections involving derivatives 
of theta functions [Eyn09]. Taking this series as a definition of an object T[^,u]{S), 
some of its properties related to modularity and holomorphic anomaly equations were 
studied in [EM08]. In this article, we argue that 7[^,i/](iS) is a good candidate for a Tau 
function. 
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5.2.1 Preliminaries 

We shall consider the dependence of the F^'s in the filling fractions^: 

= Pgi^i^ • • • ^eg) 

and we shall write for short F' for the vector of first derivatives of Fg with respect to 



e,: s: 



F' 



dFg dFg 



and Fg, Fg", . . . , Fg'^\ denote likewise the tensor of higher derivatives. 

We need Theta functions, and we shall introduce an appropriate notation for our 
purposes. We define the Q[fiM] function as: 

'^[iiu]{^\'T) = ^i^T{p+^^-Ne)■T■[p+^l-Ne)+{p+^l-Ne)■vf+2i■Kp■v j^g^ 

[/i, v] play the role of a characteristics, although we do not require it to be half-integer 
here. We can compare to the definition of the Siegel theta function with characteristics 

^[l^](w|r) = ^ exp[z7r(p + /i)r(p + yu) + 2z7r(p + /i)(w + u)] (5.17) 

and thus: 

e[^,,](w|r) = e'"^'"""-^""'-2*"^-'^i9[^;](w/2m - Nr ■ e|r) (5.18) 

We shall omit to write the dependance in [yU, v\ , and we use the notation 9', 9", . . . , 9^*^) 
for the tensor of derivatives with respect to w. 

9'= fl^,...,-^^ andsoon (5.19) 



^dwi ' ' dw 
This 'd[^,u] function satisfies the heat equation: 

1 (a.,,^, + 9.^,,, ) 9[,,,] (w) = Z7r9"(w),, (5.20) 

Again, in this equation, Tj^ and considered independent. 

From now on, the dependence in the characteristics will be implicit. Given a spectral 
curve (C, X, Y) of genus q, and the choice of a symplectic basis {A, B), we set: 

e = — (p YdX 

Wo = NF^ = N (f YdX 
Jb 

T = F!! = <j) d) B the Riemann matrix of periods of C 



^We warn the reader not to confuse the the genus g of the spectral curve, and the index g of the 
symplectic invariant Fg. 
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and we define to shorten formulas: 



0^'^ = (w = Wo I r) (5.21) 

as tfie tensor of k^^ derivatives of 0. In the definition of 0, notice that wq appears 
through the combination: 



( = Wo - iVr ■ 2nie = N (k YdX (5.22) 

Jb-tA 

which is modular of weight 1, and already appeared in Section 3.1. 
5.2.2 Definition and comments 

Definition 5.3 For any (//, u) G and any spectral curve S = (C, X, Y) of genus g, 
we define: 



+o(l/iV2)j (5.23) 



1112 




This function is defined formally order by order in N (in the coefficient ofl/N^, each 
term Q^^^ is considered formally of order 1). It can he seen as a genuine asymptotic 
series when the reference filling fractions e are chosen such that wq = NFq is of order 
0(1). 

When /i + ri/ is a non-singular half-integer odd characteristics, the leading term when 
— oo of T{S) coincides with the classical Tau function of Section 4 computed for 
the differential form u = N YdX, up to an exponential factor: 

V,](5) ~ e^'^»e(wo|r) 

~ e^'^»^(C + z/ + /i-r)e2^"^< 
~ rci{NYdX)e^'''^-^ 

Notice that the Fay identity for 6 [Fay 70] presented in § 2.2.3 is also true if we multiply 
6 by an exponential factor of its argument. Therefore, Theorem 4.3 ensures that the 
large limit of T{S) satisfies the Hirota bilinear equation (Theorem. 4.3), i.e. is a 
dispersionless Tau function to leading order. This is also true for arbitrary character- 
istics although the times have to be shifted by a (maybe complex) constant in 
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this case. We conjecture (see Section 7) that TJ^.j^] is actually a Tau function, satisfying 
Hirota equations to all orders. 

We also mention that under a modular transformation (i.e. a change of choice for 
the A and B cycles), T[fj,^u] changes like the Theta function of characteristics [fi, u] (see 
[EM08]): 

Proposition 5.1 Under a modular Sp{2Q, Z) transformation r — )■ f = {At + B){Ct + 
D)^^ , the characteristics [/i, u] changes as n ^ jl = Dfi — Cv + |(Ci5*)diag; v ^ u = 
—Bji + Au + |(y4i?*)diag, ihe Tau function transforms as 

Ti^M ^ C[m,h(^' B, C, D) (5.24) 

where ([^ ,^]{A, B,C, D) is the phase factor, independent of the spectral curve, which 
already appears in the modular transformation of the Theta function. 



5.3 Baker- Akhiezer spinor kernel 

We now define the spinor kernel ipi^iy ^2 ', 5) through Sato's relation: 
Definition 5.4 

U T{S,YdX^YdX + j^dS.,,.,) 
ij{zi,Z2; S) = ^3^^^ ^ydX{zl)dX{z2) 



nS+[z,]-[z2]) 



y/dX{zi)dX{z2) 



ns) 

where [z] is Sato's notation, see Eqn. 4-9 in Section 4-1- 

ip{zi,Z2) is again defined formally, order by order in or as an asymptotic series 

for appropriate e. The leading order coincides with ipci introduced in Definition 3.1. 
Let us give the first few orders: 

_ ^'^^ e(wo + 2iTr{u{zi) - u{z2)) \ r) 

'ip[Zi,Z2 ; oj - 



E{zi,Z2) 6(wo|r 
1 



1 



6 



(0) _L / , ,(1) 



UJ3 + UJi 

r) fzx /-zi 



+- 



1 e'(wo + 2i7r(u(^i) - \x{z2)) 



2 e(wo + 227r(u(2i) - \i{z2)) I r) 
i e^Xwo + 2z7r(u(zi)-u(z2))|r 
+ 2 e(wo + 2z7r(u(zi)-u(z2))|r ' ^ ^ 



(0) 
3 



2 J Z2 JB 

z\ 



Z2 JbJb 



e^(wo + 227r(u(2;i)-u(£2))|r) _ 9^(wo | r) , ^, 
e(wo + 2m(u(2i)-u(^2))|r) e(wo|r) ' ^ 
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0'"(wo + 2z7r(u(zi: 



6 I 6(wo + 2i7r{u{zi] 



■ "(^2)) k) 


0'"(wo 


\r) 


u(;^2))|r) 


e(wo 


r) 



^0 



+ o{l/N) 
(5.25) 



Lemma 5.1 il){zi,Z2:, S) is a well-defined spinor in Z\ and z^. Furthermore: 
• "^{zx, Z2) has a simple pole at zi = Z2: 

^dX{zi)dX{z2) 



2^1^22 X{Zi) — X{Z2) 



N ^=1 YdX 

o zn 



• It has an essential singularity near any pole p of YdX , of the form e •'^2 

• At all orders in 1/N (except at leading order), iIj{zi,Z2) has poles at the ramifi- 
cation points tti. Their order increase with the order of 1/N . 

Proof. When Zi goes around an ^-cycle, dSz^^z2 is unchanged. When zi goes around 
a cycle Bj, dSz^^z2 is shifted by a holomorphic form dSz^^z2 ~^ ^5*2^^22 + 2i7rdnj, which is 
dual to a d/dej, and since the Tau function is background independent (it was proved 
in [EM08] that dT/dei = 0), then it is unchanged. This shows that iIj{zi,Z2 ; 5) is a 
well-defined spinor. 

Then we compute each term of T{S + [2:1] — [^2]) by writing the Taylor expansion: 



Fg ( iS, ydx — !■ ydx 



E 

n>0 



(5.26) 



A=0 



which we need to evaluate at A = 1. The n-th derivatives of at A = are computed 
by the special geometry relations theorem 5.2, using the dual cycle (dS^-^^^J* = [z2, zi]: 



d^F„ 









A=0 


J Z2 


J Z2 



(9) 



(5.27) 



All the (jj^f ^ with 2 — 2(7 — n < are meromorphic, and have poles only at ramification 
points, without residues. This implies that their contribution to '^/'(-Zi, ^2) provides only 
poles at ramification points. The only terms involving uj^n^ with 2 — 25? — n > 0, are 
9aFo, dlF^ and dy^F^. 



dX 

d^Fo 



dX 



= j^^ FdX, which contributes to i}) as the essential singularity e^-^^a ^"^^ ^ 
— In {E{zi, Z2Y dX{zi) dX{z2)) which contributes to ip as 1/ E{zi, Z2). 
2?7r {vi{zi) — vl{z2)) which does not yield any singularity. 
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A=0 



A=0 



All the other terms have 2 — 2g — n < 0, and contribute order by order, only as 
combinations of meromorphic forms and derivatives of Theta functions, having poles 
at ramification points conveyed by the oon^^s with 2 — 2g — n < □ 

Notice that as a corollary of Proposition 5.1, ip has nice modular properties: 

Corollary 5.1 Under a modular Sp(2g,'Z) transformation r — )■ f = {At + B){Ct + 
D)~^, the characteristics [fi, u] changes as fi ^ fl = Djj, — Cu + |(Ci5*)diag; y ^ u = 
—Bji + Av + |(Ai?*)diag, ih^ spinor kernel ip[fj_^^]{zi, Z2) transforms as 

i'itiMi^^^ ^'i) ^ V^[A,i>](^i)^2) (5.28) 



6 Correlators 

6.1 Second kind deformations of S 

Let us recall the definition of the insertion operator 5^, already encountered in Sec- 
tion 4.3 and adapted now for varying spectral curves. 

Definition 6.1 We define the insertion operator 5z, acting on a functional f{S) of a 
spectral curve S, as follows: 

6 J = dX{z) ^ f{S,) 

A=0 

where Sx is the family of spectral curves obtained from S by: 

Biz.-) 



{YdX)x = YdX + X 



dX{z) 
6z is a derivation. 

The dual cycle of B{z, ■)/dX{z) is the contour surrounding z with index 1: 

Then, the relations of special geometry (theorem 5.2) for Un^ implies, for any n,g > 0: 

SzUJ^f^zi, . . . , Zn) = Res Unliiz' , Zi, Zr, 



Ujl^L{z,Zi,...,Zn) (6.2) 



{X{z')-X{z)) 

For instance: 



6zFo = u'^'iz) = Y{z) dX{z) , 6zFg = u['\z) , M '(^') = B{z, z') (6.3) 



37 



6.2 Definition 

Definition 6.2 For n positive integer, we define the correlators Wn{zi, . . . , Zn) and 
the disconnected correlators Wnizi, 



Wn{zi, ...,Zn) 
Wn{Zl, ...,Zn) 



N-^ 5,, ■ ■ ■ In nS) 



Wn{zi, . . . , Zn) and Wn{zi, . . . , Zn) are (1, . . . , l)-forms, symmetric in their n variables. 

n times 

Each coefficient, order by order in 1/A^, is a meromorphic form with poles at ramifica- 
tion points 



6.3 Examples 

For instance the three first orders of Wi are: 
Wi{z) = NY{z)dX{z) + {lney -linduiz) 



-0^ 



+o(l/iV) 
where we recall that: 



z ' j,k 



S,FS _ 1 

2i7T 2iTi 



(0), 



AzTT V Res K{z,z')duj{z')duk{z') 

^^'^a, dX{z')dY{z') 
y Res K{z,z')B{z' ,'^) 



227rdu(z) J 
(6.4) 



z'-^ai 



(6.5) 
(6.6) 



For the 2-point correlator, the three first orders are: 
W2{zi,Z2) = 5(zi,Z2) + (lnG)"- 2i^du(zi)®2i7rdu(z2) 



+ 



+ 



(i^i)'^ 



0"w 



■zi, -^2) + ( -Q- ) • [^7r(5^iT) ® 2i7rdu(22) + 2iTtd\i{zi) ® ^-^-(J^jr)] 



e G2; ^~6"v~e~ 



0'"0'w 



02 



2i7rdu(zi) (8) 2i7rdu(z2^ 



+o(l/iV) 



(6.7) 
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For n > 3, the leading order of the ra-point correlator is a 0(1), and is obtained by 
successive applications of 6z to the (In 0)" term. 



Wnizi,--- ,Zn) 



(lne)(") •(g)2i7rdu(zj) 

'e"\("-i) 



+ Yl (lne)("^i)-( (•,z,-,Zfc))»(g)2i7rdu(zO 



[iirdzj) ® 2i7r du(zfc) 



l<j<fc<n 



+ 



6 V e e2 



2i7r du(zj)| 



-o{l/N) 



(6.8) 



6.4 Loop equations 

Theorem 6.1 The dispersive Tau function obeys the loop equations. Namely, let us 
denote hyV a contour separating the poles ofYdX, from the set of preimage of a point 
xeC. 

6zT[S]-N{YdX){z) 



zer X{z)-x 







/ A (^z6z'T[S] - dX(^)dX(^0 \ ^ g^^) 

Ler^'^^ {X{z)-x){X{z')-X{z))\ {X{z) - X{z')y J 

where Q is a rational function ofx, whose only poles are those ofYdX , and with degree 
one less than that ofYdX. 

Those loop equations can be written in terms of correlators by applying ' ^z„- 
Theorem 6.2 Let J = {z2, . . . , Zn}- The correlators satisfy 
(i) The linear loop equations. For all n > 1: 

L xw^ ("'"'^' - " - wl0Mr) ^ ° 

(ii) The quadratic loop equations. For all n> 1: 

± w„,,(., J) + dx(.) d^(.') 5: d,. ( (^_;^;'//j^(,.^) ) } = J) 



39 



defines a quantity Qn{x', J) which is a rational function of x, whose only poles 
are located at those ofYdX. 

The important information in loop equations, is that those particular combinations of 
WnS have no monodromies in the variable x. Here in particular, they are regular at 
the ramification points aj. Since every Wn has poles at ramification points to all orders 
in this is a highly non-trivial property. 

Proof. The formula: 

T{S) = J2 e"'^""'" Z[S,Yi] (6.9) 

makes sense, provided we collect all terms of a given order in — )■ oo as in Defini- 
tion 5.3. The Fg were precisely introduced such that for any n, 

Z[5,n] =e^«^'"''^«('5'"+^) (6.10) 

is a solution of the loop equations. By linearity, T{S) satisfy the same loop equations. 
In fact, the definition 5.3 of T{S) was introduced in [Eyn09] precisely to give a solution 
to those loop equations, which has the extra property of being independent of a choice 
of filling fractions. □ 



7 Hirota equations 

In the dispersionless case, we have seen that the self-replication property for ip^x is 
equivalent to an infinitesimal Fay identity, which is known in turn to be equivalent to 
Hirota equations. In a similar way, we conjecture here: 

Conjecture 7.1 ip is self-replicating: 

(5^^)(2i, Z2) = -ipizi, z)i){z, Z2) 

We have not been able to prove this conjecture. We prove in appendix A that it holds 
up to o(l/A^). We argue in Section 9 that it is compatible with what is known for 
spectral curves coming from the one matrix model (hyperelliptic curves), or the two 
matrix model. Besides, these matrix models do not allow to reach all plane curves 
S. The difficulty in finding a proof of Conjecture 7.1 comes from the singularities at 
ramification points. For instance, one can always write: 

4 {6,ij){z^, Z2]S)= Res ij{z, z' ; S) ^{z' , ^ ; 5 + [2^1] - [2^2]) (7.1) 
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where S + [zi] — [Z2] is the spectral curve S where YdX is replaced by YdX + jj dS'^^^^j- 
Since the integrand is a differential form on the Riemann surface S underlying S, we 
can move the contour to the poles at the ramification points, and the pole at z' = Z2: 

^ {6,ij){zi,Z2; S) 

= -^{zi, z ; S)tp{z, Z2; S) -'Y] Res tp{z, z' ; S) tp{z', z; S + [zi] - [Z2]) (7.2) 

i 

Then, it remains to show that the sum of residues at ramification points vanishes. So, 
conjecture 7.1 is equivalent to: 

Conjecture 7.2 

V Res iIj{z, z' ■ S) ijiz', z;S+ [z^] - [Z2]) = 0. (7.3) 

i 

In Appendix A, we check that this residue at each is o(l/A^). This involves already 
non-trivial identities between Theta functions associated to a Riemann surface, like Fay 
identity or its degenerations, and involves the precise value of uj^^ ■ We have yet not 
been able to find a general way to show that this residue is to all orders in 1/A^. 

In terms of the T function. Conjecture 7.1 can be rephrased: 

Conjecture 7.3 T satisfy a "differential Fay identity": 

r[s] {5,r)[s+[z,]-[z2]]-{8zr)[s] r[s+[z,]-[z2\] = -Nr[s+[zr]-[z]] r[s+[z]-[z2]] 

There is also a global version of the former conjecture. First, notice from our definition 
in Section 5 that: 

T[{S + [Z,] - [Z2]) + [Z,] - [Z,]] = r[iS + [Zs] - [Z,]) + [Z,] - [Z2]] 

so omitting the parentheses makes sense, but: 

T[S + [z,] - [z,] + [^3] - [^2]] = -r[S + [z,] - [Z2] + [z,] - [z,]] 

This sign comes from the fact that the definition of T contains a regularization pro- 
cedure (for f^^ f^^ B), whose result depends on the way we form the pairs of simple 
poles to add to S. 

Conjecture 7.4 T satisfies a "Fay identity": 
T[{S+[z,]-[z2]) + [z,]-[z,]]T[S] 

= r[s + [z,] - [Z2]] r[s + [zs] - [z,]] - r[s + [^3] - [^2]] r[s + [z,] - [z,]] 
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Provided our conjectures hold, T and are actually the Tau function and the Baker- 
Akhiezer spinor kernel of a dispersive integrable system. 

Proof of equivalence of Conjectures 7.3 and 7.4. 

We can obtain Conjecture 7.3 from Conjecture 7.4 by letting zi and Z2 merge to a point 
z. In the other direction, we use condensed notations 

%jki = r[s + [zi] - [z,] + [zk] - [zi]] , 71, = r[s + [z,] - [z^]] , r = r[s] 

and we apply Conjecture 7.3 to the spectral curve S + [z^] — [z^]: 

Oz 11234 — — /l234 ~^V — 

/34 / 34 

^zT ^ ^jT3zTz4Tl234 ^jTlz34Tz234 ,-. .n 

- /1234 - — — iV — {(A) 

I I 134 134 

Exchanging the roles of Zi -H- 23 and Z2 -H- 2:4 also gives 

r ^ ^zT ^ ^jTlzTz2Tl234 j.rTszu'T'ziU ,^ r\ 

0zTl234 = '1234 " N — — — (7.5) 

/ / 112 112 

and comparing the two, we may get rid of the terms involving 5z'- 

Tl273zTz47l234 + '7~7l27l234 7^234 = 'T2,a'TizTz2'Ti234 + TT34'T-3zl2Tz412 

(7.6) 

Let us define 

f^l234 = 7~Tl234 — '7l2'734 + TliTi2 ■'^) 

which is the quantity which should vanish at the end of our computation. Notice that 

lim^^^22Wi234 = (7.8) 

Indeed: 

^1234 ~zi^Z2 ~^7~ h T5l734 ^7- 5iT T-3A + TliT3l 

J^12 J^12 

^zi^z2 T61T34 - S1TT3A + T14T31 (7.9) 

where Eij = E{zi,Zj), and this expression vanish by application of Conjecture 7.3 to 
the spectral curve S. Notice also that the remark about T made above Conjecture 7.4 
implies U^ki = Ukuj = —l^ukj- Let us rewrite Eqn. 7.6 in terms of Uijki only: 

712^1234^2234 — 734^3212^2412 + {TuTszTzA — TlzTz2T3A)hli23A 
+ {TlzT34 — Tl4T3z)Tl2Uz234 + (7^2734 — Ti2TzA)Tl2Ulz34 
+ (732 7l2 — Ti2T3z)T3aUz412 + (7^27l4 — 7127^4)734^3212 
= 
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The left hand side may have simple poles when Zi — >■ Zj due to 7ij, but not higher degree 
poles since we know that lAijki is actually regular when Zi — )■ Zj or zi. For Eqn. 7.10 to 
hold, in particular, the coefficient of the simple pole when Z2 — )■ z^ must vanish: 

- Tn%A + Tu%A W.413 = (7. 10) 

which we can also write after reindexing the points and using the symmetries of U: 

TziTzaUizZA = TzzTiaUzIM (7-11) 

Similarly, the coefficient of the simple pole when zi — )■ Z2 must vanish: 

^1^34^2134 + (7^l734 — 7^473l)Wi234 + (71^734 — 714732)^2134 = 

(7.12) 

Now, we may combine the latter with Eqn. 7.11, set 2; = 22 for convenience, and isolate 
W1234: 

-7- -7- ^^1234 
/32 /14 

((r2ir34 - r24r3i) + ^^^{ri2%A - riA%2))u^2ZA = o (7.13) 

^ /32 / 14 ^ 

If W1234 were not identically zero, we would have (by continuity of all the coefficients 
of the series) for any points Zi,Z2,z^, z^ on the curve: 

732 7l4 

W1234 = _ _ (724731 — 721734) + (7i4732 — 712734) (7-14) 

/3I '24 

But matching the coefficient of the simple pole when zi — )■ 2:4 in this equation yields: 

= -^(r2ir3i - r2ir3i) + r32 = r32 (7.15) 
/31 /21 

which is not true. Therefore, W1234 = 0. □ 

8 Consequences 

8.1 Exponential formula 

We start with a remark which does not rely of the conjectures of Section 7, but which 
is natural to present now. Recall that the kernel is defined by Sato's formula: 

II r(5,ydx^ydx + ^dg2,..J . 

ij{zi, Z2;S) = ^3^^^ '-^ V^^i^i) dX(2;2) (8.1) 
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Adding a double pole to YdX can be realized by adding two simple poles and take 
the limit where the two simple poles collapse. In other words, we can write dSz^.zj = 
= {X{zi)-X{z2))B{.,zi)/dX{zi) + 0{{X{zi)-X{z2)))\ and then express ^ in 
terms of second kind deformations of T, i.e. in terms of the correlators Wn- However, 
when we substitute second kind deformations instead of third kind deformations, we 
must pay attention to regularization for the term B, like in Section 4.1 and 5.3. 
Besides, for any n, Wn has a finite number of terms which are not 0{1/N). For this 
reason we define the following quantities, which are 0(1/ N): 

Wi{z) = Wi{z) - NYdX{z) - {\ney ■2i7rdu{z) (8.2) 
W2{zuZ2) = W2{zi,Z2)-B{zi,Z2)-{\ne)" ■2i7rdu{zi)®2iTTdu{z2) (8.3) 

and for n > 3: 

n 

Wn{zi, ...,Zn) = Wn{zi, . . . , z„) " (In e)^") ■ (g) 2z'K du{z,) (8.4) 

Then, we have the exponential formula: 
Proposition 8.1 

E{zi,Z2) e V^n!^ ^ 
where 

812 = e(wo + 2i7r(u(2i) -u(^2)) |r) , e = e(wo|r). (8.5) 
This formula is an equality if we collect on both sides all the terms of the same order. 

Proof. The Taylor formula allows to express \nT[S + [zi] — [Z2]] — lnT[5] to all orders 
in For the first and second order, we use Eqn. 4.2 and Eqn. 4.3. 

lnT[S+[zi]-[z2]]-\nT[S] 

"21 1 pZl pZl 



r N-H^r^T) + \ " N-\6^.{\nT)[S] " 

J Z2 J Z2 J Z2 

M~n rzi pzi 

„^n J ZO J ZO 



' Z2 



n times 
21 /"Zl 




n 

n>3 
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N ^ ' YdX - i In ({E{z,, Z2)fdX{z^)dX{z2] 

n>l 

+ 5:-/ .../ (8.6) 



77 ' 

n>l C^^2 "^^2 
n times 



In the last step, we have used the expression for the leading order of Wn found in 
Section 6.3. Then, the second line can be resummated into In G12, and the whole result 
exponentiated leads to the formula we announced. □ 

Notice that, when expanding the exponential, to any given order 0{N~^), all the 
Wn give a contribution involving derivatives of theta functions contracted with tensor 
products {u.(zi) — u(z2))^"'. These contributions have to be resummated into a single 
theta function (or derivatives thereof) with argument shifted by u(zi) — u{z2), then 
producing an expression at the order sought involving only a finite number of terms. 

8.2 Determinantal formulas 

Conversely, the correlators Wn{zi, . . . , Zn) can be recovered from the spinor kernel 
ip{zi,Z2), they are the determinantal correlation functions built with ip{zi,Zj). Let us 
consider first the case of Wi, which does not rely on the conjectures of Section 7. 

Lemma 8.1 

W,iz) = NYdXiz) + hrn^\^i;{z,z)e " V) - X(.) j 

Proof. First, we notice that adding a double pole to YdX can be realized by adding 
two simple poles and take the limit where the two simple poles collapse. More precisely, 
when z' ^ z, we have 

dS,,,{zo) ~ (X(/) - X{z)) (8.7) 
z'-^z aX [z) 

We can thus use the definition of 6z with A = X{z') —X{z). For any regular functional 
f(YdX) of the spectral curve, we thus have: 

/(ydX + ldg...) = fiY^X) + ^^^\x{z)'^ 6J{YdX) + 0{{X{z')-X{z)r) (8.8) 

In some sense we trade a variation of YdX by a second kind differential with a variation 
with a third kind differential. 
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In particular, Fg with > 1, or every 0-term in Definition 5.3 are regular func- 
tionals of the spectral curve. We just have to pay attention to the Fq term, because 
the derivative of Fq with respect to third kind differentials involves a regularization 
procedure, whereas the derivative with respect to second kind differentials does not. 
We have, by Taylor expansion, and computing all derivatives from special geometry 
(Theorem 5.2): 



Fo(rdX + Ad5, 



Fo(FdX) + FoiYdX + A dS, 

n=l 

Fo{YdX) + xj^ YdX 
-yln {E{z,z'fdX{z)dX{z')) 



n=3 

Fo(FdX) + A I YdX 

J z 

-yln {E{z,z'fdX{z)dX{z')) 
+0{{X{z') - X{z)f) 



A=0 



Taking A = gives: 



1 



FoiYdX + - dS,, + iE{z', z) VdX(z)dX(z')) 

+0{{X{z')-X{z)r) 



1.10) 



Finally we have: 

tljiz\ z) e-^/^^' E{z', z) = l + ^^t}'^^''^ 5, (In T - N^F^) + 0{{X{z') - X{z)f) 



NdX{z) 



(8.11) 



and thus 

i'{z\z)e 



£/(r,2;) A/ dA (z) ii/(2;', z) 

(8.12) 



Taking the limit z' -> z, and noticing that l/E{z',z) = y/dX{z) dX{z')/{X{z') 
X{z)) + 0{X{z') - X{z)), gives the lemma: 

1 



6,{\nr-N'Fo) = Wi{z) - NY{z)dX{z) 
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□ 

Theorem 8.1 If Conjecture 7.1 holds, then 

n 

\/n>2, Wn{Zi,...,Zn) = {-ir+' 

(T cyclic perm. i=l 

Equivalently: 

Wn{zi, . . . , 2;„) = " det " tlj{zi, Zj) 

where " det " means that when we decompose the determinant as a sum of permutations, 
each factor tp{zi, Zi) should he replaced by Wi{zi). 

Proof. The formula for Wi is proved in Lemma 8.1. Then, we get the formula for Wn 
by recursively applying 5^- and using the self- replication of ^Z^. □ 

It is thus clear that this determinantal structure relies on the (conjectured) existence 
of Hirota equations for T, i.e. on integrability. 

8.3 Baker- Akhiezer functions 

We recall the notations of Section 3. With the kernel iIj{zi,Z2) of our dispersive con- 
struction, we introduce a. d x d matrix \Ef(xi,X2 ; S) = ipiz'^i^i), z^ {X2) ; S)ij=i^,,,^d for 
xi,X2 G C This definition has to be regularized when xi or X2 is equal to X{p) where 
p eV. In particular, it was explained in Section 3.4 how to take xi or X2 — )■ 00: 

^{x-S) "=" [^(2^(x),oo,;5)],,i<,<rf (8.14) 
$(a:; S) " = " [^(cx),, z^(x)) ; S]j,i<,<d (8.15) 

8.4 Duality equation 

Theorem 8.2 // Conjecture 7.1 holds, we have for any spectral curve S: 
^(xi, X2 ; S)^{X2, xs;S)= ~ -^g) dxa ^^^^^ _ 

[X2 — Xi)[X3 — X\) 

Proof. For S = {C,x,y), we introduce an auxiliary spectral curve: 

4- =(c,x,Y + ^ ^^f:M^\ (8_ig) 
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To compute the matrix element, we use the self-rephcation of ijj: 



Now, we use the fact that 6z is a derivation, and by definition, 6z lnT[4] — ! <Sij]- 
[^(xi,X2;5)*(x2,X3;5)],, = - ^ M/i[^™(x2) ; 4] V^(^'(xi), ^^■(xa) ; 5) 

m. 

+ J2Wi[z'^{x2);S]ij{z\x,),z^{x3)-S) (8.18) 

The hnear loop equation (Theorem 6.1) tells us the sum over sheets of Wi[z"^{x2) ; 4]- 
[^'(xi,X2; 5)^'(x2,X3; S)]ij 

= [Vl/(X1,X3; 5)],,- (- ^iVFdX(z-(x2)) -d5,,(,,),,.(,3)(z'"(x2)) 



+ A^FdX(2"^(x2)) 
= -[^(xi,X3; S)]ij ^dS'^i(^^),^,(^3)(z'"(x2)) 

m 

(Xi — X3) dX2 / / N 

(X1-X2)(X2-X3) 

□ The spinor kernel ipci{zi, 22) of the 
dispersionless integrable system of Section 3 was regular at ramification points. So, we 
could find a formula for ijjf.i{zi, z)iIjc\{z, Z2) even before the sum over the sheets where 
z is located (refined duality equation. Theorem 3.3). Here, the spinor kernel iIj{zi,Z2) 
does have, order by order in 1/A^, poles at ramification points. So, we do not have a 
simple expression for ip{zi, z)ip{z, Z2)- However, Theorem 8.2 shows all contributions 
from the ramification points cancel in the sum over sheets. 

8.5 ChristofFel-Darboux relations 

The matrix \l/(xi, X2 ; S) is invertible, since it is defined as a series for which the leading 
term coincides with \l/ci(xi,X2; S) which is invertible (see Lemma 3.2). We also have 
a duality relation to express the inverse: 

Corollary 8.1 If Conjecture 7.1 holds, 

^-\x^,X2 ; S) = ^(X2, xi ; S) 

dxidx2 
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Theorem 8.3 // Conjecture 7.1 holds, the matrices S) and S) are invert- 

ible, and we have the Christoffel-Darhoux relation: 

• = xW^) 

where the matrix A is invertible, independent of x and given by: 

A-^\S] = S)^'{x; S) 

ax 

The proofs can be copied from the classical case (Section 3.5), because they are based 
only on the duality equation. 

8.6 Differential systems 

\&(xi,X2; S) is the solution of a system of differential equations with respect to the 
positions of the poles X{p) (x = oo is a fixed pole) and the times tpj. 

Theorem 8.4 // Conjecture 7.1 holds, then for any deformation parameter A = tpj 
or X = X{p) 7^ oo, there exists a d x d matrix Mx{xi,X2 ; S), such that: 

(i) (9a-Ma)*(xi,X2; 5) = 0. 
(ii) M\ is a rational function of X2 ■ 
{Hi) Mx has no pole at branchpoints. 

These deformations are compatible (since \I/ is invertible) and isomonodromic. 

Implicitly, for j = 0, only deformations dx = ^ — (9y ^ are considered. 
Proof. By definition: 

Mx = dx^{Xi,X2)^^\xi,X2) 

= - ^'''~r^' dx^{x„X2)^{x2,x,) (8.20) 
When A is a time tp^k {k > 1), we have by self-replication: 



dt^,M^,,Z2) = - I iJ{z^,■)H■,Z2) (8.21) 

J* , 

p,k 



We compute: 



1 f {X{-) - Xi){x2 - Xi] 



da^i L;, (^(0-^2) 
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■ij{z\,-)ij{;zi) (8.22) 



It is clear that Mt^^,{xi, X2) is a rational function of X2 without poles at branchpoints. 
If we assume X{p) 7^ 00, it has a pole only at X2 = X{p), and this pole is of order 
k + 1. If we assume X{p) = 00, it has a pole only at X2 = 00, and this pole is of 
order 1 + [k/dpl (we recall that dp is the multiplicity of order of the pole p of X). It 
is straightforward to adapt this proof for k = 0, namely dx = dt^Q — dt^, ^. 

Now, we turn to \ = X{j>) ^ 00. We have assumed that dX{p) 7^ 0, so the 
preimages p^ , . . . ^p'^ oi X{p) are distinct, and {X{z) — X{p)) is a local coordinate near 
each p™. The Laurent expansion of YdX{z) when z ^ p^ is: 



y<m^)-l,t,. ^^J%^Y,. +0(l) (8.23) 

where only a finite number of tp^k' = are non zero. If we perform X{p) X{p) + X 
while the times are fixed, we change YdX to (YdX)x, with Laurent expansion at p"*: 

iVaX), = ydX + A ^ + Oil) (8.24) 

So, we can identify: 



d 



dX{p) 

which is a finite sum. Hence: 




(8.25) 



k>l 



This is a rational fraction of X2, with a pole at X2 = X{p) of order 1 + 
[(maXmOrd2m(p)FdX) /dp\. □ 
One could also define a matrix L with the reconstruction formula presented in 
Section 3.6. But now, the function Y depends on the spectral curve, the evolution of L 
under the flows generated by the times is not isospectral. But the argument of Eqn. 3.8 
is still valid, and the times in this construction are isomonodromic parameters. Such 
a reconstruction has also been performed at leading order by Bertola and Gekhtman 
[BG07], and our expressions up to o(l) match indeed their results. 



9 Parallel with matrix models 

Some unitarily invariant models of random matrices are special cases of this construc- 
tion, that actually suggested to develop this formalism. We recall here the dictionary 
between the definitions in integrability (and in this article), and observables in matrix 
models. Let M be a random square matrix of size N x N, diagonalizable by a unitary 
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conjugation, with eigenvalues restricted to some contour T in the complex plane. The 
probability measure d/i(M) is a data of the matrix model. Its normalization defines the 
partition function Z^r = / d/i(M). For a function /, we note (^/(M)^ the expectation 
value of the random variable /(M): 

9.1 Examples of integrable matrix models 

Important examples of integrable matrix models, described in details in [Meh04], in- 
clude: 



The one hermitian matrix model 

M G H]\f is a. N X N hermitian matrix, with eigenvalues on the real axis 

N 

d/i(M) =e-^^^(*^MM , dM = Y[dMi,i Yl dRe Mi^j dim Mi J (9.2) 

i=l l<i<j<N 

and V{M) is a " semiclassical potential" (see [Ber07]), i.e. V'{x) G C(x), i.e. V'{x) 
is a rational function of its variable a;, chosen such that e~^^^^dx is absolutely 
convergent. 

The one normal matrix model, with eigenvalues on a contour F 

We define 

Hn{T) = {Udia.g{xi,...,XN)U^ t/ G t/(iV) and (xi,...,X7v) e F^} (9.3) 
equipped with a " measure" (not necessarily real or normalized) ' : 

N 

dM = dU Y[{xi - Xjf Yl dxi (9.4) 

i<j i=l 

where dU is the Haar measure on the unitary group U (N) and dxi is the curvilinear*^ 
measure along F. This measure is always invariant under unitary transformations. We 
then define the measure d/i(M) on ifjv(F): 

d^M) =e-^'^^(*^MM, (9.5) 



^For a given M S HN{r), there exist many choices of U and z^'s (in fact U can be multiphcd by 
any element of U{1)^ and the Xi's can be permuted, but the reader can check that the measure dM 
is well defined on Hn{T). 

^If r C C is a path in the complex plane parametrized by a function 7 : M F C C, 
i.e. F = {7(5) , s € M}, at X = 7(5) we define the curvilinear measure dx ~ j'{s) ds where ds is the 
Lebesgue measure on M. 
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where V'{x) G C(x), i.e. V'{x) is a rational function of its variable x, chosen such that 
Jp e~^^'^-*dx is absolutely convergent. This means that for a given potential V{x), T 
must go to oo or to the poles of V, only in sectors where KeV ^ +00, i.e. in some 
"Stokes sectors". When F = M, this definition is the same as the Hermitian matrix 
model Hi^[M.) = Hf^r, with the usual Lebesgue measure on if^r. When F is the unit 
circle Si, this definition is the same as the unitary matrix model Hn{Eii) = U{N), with 
the usual Haar measure on U{N). 

The two normal matrices model 

Given two paths Fi, F2, we define a measure on i77v(ri) x -f^Ar(r2) 

d/i2(Mi, M2) = e-^T-- mM,)+V,(M,)-M,M,] ^9_Q^ 

where Vi and V2 are semiclassical potentials {V( and V2 are rational functions), chosen 
such that integrals on Fi x F2 are absolutely convergent. Upon integration on M2, this 
measure induces a measure dyu(M) on Mi that we rename M: 

d/i(M)= [ d/i2(M,M2) (9.7) 

The chain of matrices 

This is the natural generalization of the case of two matrices. Consider k paths 
Fi,...,Ffc, and k semiclassical potentials Vi,...,Vfc, and define the measure on 
HNiTi) X ... X HNiTk) as 

d/i,(Mi, . . . , Mfc) = e-^^^ E'- y^iM.)-j:^Zl A/.A/.+i] (9.8) 

For any i G {1, . . . ,k}, we integrate on Mj's with j 7^ i, and renaming M = Mi, this 
measure induces a measure d/i(M): 

d/i(M,) = / d/ifc(Mi,...,Mfc). (9.9) 

9.2 Correspondences 
9.2.1 Partition function 

All the listed matrix models above have the property to be integrable, in a sense 
explained below. There is a huge literature on the subject, let us mention among 
others the early works in physics [MMKZ92, Kos96], and in mathematics [HTW93, 
AM95, vM97]. For the one and two matrices models, it is established that the partition 
function Zn = J d/i(M) is an isomonodromic Tau function [BEHOSc, BEH06, BM09a] 
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(such a result is not known at present for the chain of matrices). The so-called double 
scaling limit of matrix models has been also intensively from the point of view of 
isomonodromic deformations [IKF90, DS90, Moo90, FIKN06]. 

9.2.2 Correlators, spectral curves, loop equations 

For all the matrix models above, the correlators are: 



n _ 

W„(xi,...,a:„) = iV""(nTr--^) (9.10) 

i=i * 

n _ 

Wn{xi,...,Xn) = N-^'^HTi^j^) (9.11) 

1=1 "^^ 

where C means 'cumulant'. Wn is often called n-point function, or connected n-point 
function, and W„ is called disconnected n-point function. 

The one point function >Vi(x) plays an important role, it is often called the resol- 
vent. In all the matrix models above, under reasonable assumptions, it has a large 
limit denoted: 

^Wi(x) ~ Wi'\x) = y{x) (9.12) 

which is furthermore an algebraic function of x, i.e. there exists a polynomial S{x,y) 
such that 

£{x,y{x)) = 0. (9.13) 

This algebraic equation defines an algebraic curve C, i.e. a compact Riemann surface 
C and two meromorphic functions X,Y : C ^ CP^, such that 

V^gC, y{X{z)) = Y{z). (9.14) 

The " semiclassical spectral curve" S is then 

S = {C,X,Y). (9.15) 

The correlators corresponding to Definition 6.2, are 



Wn{zu...,Zn) = N^Wn{X{zi),...,X{Zn))QC)dX{Zi)+dn,2 



i=l 



{X{z,) - X{z,)Y 

(9.16) 



It is a classical result of random matrix theory (for instance it can be proved by 
integration by parts in the matrix integral) that, for all matrix models listed above, 
the correlators Wn satisfy the loop equations of Theorem 6.1. 
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9.2.3 Large asymptotic expansion 

It is conjectured that the large asymptotic expansion of the partition function 
matches with that of the T function for the semiclassical spectral curve S introduced 
in Definition 5.3: 

^iv = 7[^,.](5) (9.17) 

i.e. has an asymptotic expansion 

1/9' (p,"i\ 
In = N^Fo + Fi + In e + - f F^- + F^"— J + o{l/N) (9.18) 

The characteristics v\ of the Theta function is determined by the choice of the 
integration contour F. This conjecture was derived heuristically in [Eyn09]. 

When the semiclassical spectral curve is of genus 0, there is no Theta function and 
the expansion involves only powers of 1/A^ (in fact, powers of 1/A^^). This happens for 
the so-called "one-cut regime", and for the one matrix model, the existence of such an 
expansion has been proved for the one matrix model for real- valued, analytic potential 
[APSOl], and then the coefficients are necessarily given by the symplectic invariants 
of the semiclassical spectral curve [Eyn04]. Beyond the one-cut regime, the Riemann- 
Hilbert steepest descent analysis [IKF90, DZ95] has been applied to find explicitly the 
asymptotics up to o(l) in the one hermitian matrix model with real- valued, polynomial 
potential [BI99, DKM+99]. It features in general a pseudo-periodic behavior with A^, 
encoded in the Theta function, but the one-cut regime can also be recovered with 
this Riemann-Hilbert method [EM03]. These results have later been extended to the 
one normal matrix model with complex- valued polynomial potential [BM09b]. This 
method can be used in principle to find recursively the subleading orders, although it 
does not allow to write the answer a priori explicitly to all orders. Definition 5.3 is 
expected to give the correct answer to all orders. 

9.2.4 Baker- Akhiezer spinor kernel 

The spinor kernel is related to the expectation value of ratios of characteristic polyno- 

IXX13j1s ' 

/ det(Jf(,-i) - M) \ ^iX(z,)iX(z^) 
= ( det(XW-M) ) X(z,)^X{z,) 

The dispersionless spinor kernel ipc\{ziiZ2) is the large limit of ip{zi,Z2)- When 
sending X{z2) — J- oo, one gets, after proper renormalization that the Baker- Akhiezer 
function is 

ij{zi) = "V(^i, oo)" = (det(X(^i) - M)\ (9.20) 
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which is clearly a polynomial in X{zi) of degree N. It is a classical result of random 
matrix theory [Mch04] that this is the orthogonal polynomial of degree A^, orthogonal 
with respect to the measure dfi for matrices of size 1. The dual Baker Akhiezer-function 
for matrix models is 

^<'-' = "^<°°-'->" = ( deMA-(.VM) ) 

It can be derived from orthogonality relations of the orthogonal polynomials that those 
spinor kernels satisfy Hirota equations, also called "determinantal formula", for in- 
stance: 

(xi — X2)(x3 — X4) / det(xi — M) det(x2 — M) 



(xi — X3)(xi — X4)(x2 — X3)(x3 — X4) \det(x3 — M) det(x4 — M) 
1 /det(xi-M)\ /det(x2-M)\ 

(xi -X3)(x2 -X4) \det(x3 - M)/ \det(x4 - M) / 

1 /det(xi-M)\ /det(x2-M)\ ^2) 



(xi — X4)(x2 — X3) \det(x4 — M) I \det(x3 — M) 

In that case, is an integer (the size of the matrices), and Hirota equation is an 
equality between sequences indexed by A^, not only of formal asymptotic series. This 
relations have been proved in [FS03] for the one matrix model, and in [Ber, AP04] for 
the two matrices model, from which the case of chain of matrices can be deduced. It 
shows that, modulo the heuristic derivation of asymptotics of matrix integrals, Hirota 
equation as stated as in Conjecture 7.4) hold for the semiclassical spectral curves of 
matrix models (however, not all algebraic curves fall in this class). 

9.2.5 Differential systems 

The orthogonal polynomials, as well as their duals, in all cases above, do satisfy some 
ODE of order d [d depends on the degrees of potentials). As we have seen in Section 8, 
the proof is a consequence of Hirota equations. For instance for the one matrix model, 
this is a second order ODE [d = 1). For the two matrices model with potentials Vi 
and V2, we always have = 1 + deg (in case is a rational function, deg is the 
sum of degrees of all poles). In all cases we have a rf- dimensional vector: 

^{z) = "^^"'^ j (9.23) 

where the first entry is ip{zi). The other entries are obtained from ip{zi) with a pro- 
cedure described in [BEHOSa] and very similar to Section 3.4. This vector satisfies an 
ODE 

j^Ji^) = -^{Xiz)) ^{z) (9.24) 
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where T>{x) is a d x d matrix, whose entries are rational functions of x, that depend 
imphcitly on A^, on the coefficients of the potentials, and on the choice of integration 
contour F. The locus of eigenvalues of 'P(x), i.e. the polynomial equation 

Sm{x, y) = det {y Id^xd - T){x)) = (9.25) 

defines the finite spectral curve. The semiclassical spectral curve corresponds to its 
large limit: 

8{x,y) = lim £N{x,y). (9.26) 

There are also dxd differential systems for derivatives with respect to all coefficients of 
the potentials, and there is also a linear recursion relation on — + 1 (see [CI97] for 
the one matrix model). All these systems are compatible as shown in [BEH02] in full 
generality for the chain of matrices. Moreover, in all these matrix models, there is also 
a compatible Toda chain recursion equation relating Z^v+i, Z]\r and Zj^^i [GMM+91]. 
We postpone to a future work the investigation of the presence of such relations in our 
construction. 



9.2.6 Symplectic invar iance 

Notice that, in the 2-matrix model 

Z = [ d/i(Mi,M2), 

we have defined our semiclassical spectral curve S = (C, X, Y) from the large A^ limit 
of the resolvent VVi(x) = ( tr ) associated to the matrix Mi. Since Mi and M2 play 
a symmetric role, it is clear that we would have obtained the same partition function, 
starting from the semiclassical spectral curve S = {C, X, Y) associated to the resolvent 
of matrix M2, and thus we must have 

r{S) = r(5). (9.27) 

One can easily find that the two spectral curves S = {C, X, Y) and S = (C, X, Y) are 
related by C = C and X = Y, Y = X, in other words they are symplectically equiva- 
lent, and the fact that T{S) = T{S) can be seen as a consequence of the symplectic 
invariance of the Fg^s. 

In fact, this is a manifestation at large A^ of an exact result for finite A^. The 
orthogonal polynomials \l/(x) associated to matrix Mi satisfy an ODE of some order d 
(see Eqn. 9.23): 

-4-T^(^) =V{Xiz))^iz) (9.28) 
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whereas the orthogonal polynomials "^{y) associated to matrix M2 satisfy another ODE 
of some order d (in general d ^ d): 

j^^iz) = ViY{z)) l(z) (9.29) 

It was discovered in [BEHOSb] that: 

det {y Iddxd - V{x)) = det {x Idj, j - V{y)). (9.30) 

This also implies that the semiclassical spectral curve defined from Mi or from M2 are 
related as explained above. 

10 Conclusion 

For integrable systems with a small dispersive parameter 1/A^, using the theory of sym- 
plectic invariants [EO07a], we have introduced a formal object T, which is conjectured 
to be a Tau function in the sense that it satisfies Hirota equations. It is challenging to 
find a full proof of Conjecture 7.4. 

This construction was strongly motivated by results known for matrix models, but 
developed independently. In some sense, it gives the answer to all order in of a 
Whitham averaging. All properties mentioned in this article, correspond to well-known 
or conjectured properties of matrix models. Our construction is thus an axiomatized 
extension of the properties of matrix models, to arbitrary algebraic spectral curves. One 
can wonder how to generalize the construction. First, to algebraic curves in C* x C*, 
having in view the mirror curve appearing in Gromov-Witten theory, which all are 
of the form Pol(e^, e^) = 0. Second, to the bundles appearing in generalized matrix 
models like the 0{n) model, which may be related to isomonodromic deformations on 
Riemann surfaces of positive genus and to Hitchin systems [Hit87]. And, even further 
to P-modules, for which an adapted topological recursion (the so-called /3-deformation 
of the topological recursion) has been found with similiar properties [CEM09, GEM 10]. 

Eventually, it remains to compare this construction with other approaches (Frobe- 
nius manifolds, Poisson bracket structures, Segal- Wilson formalism in the Grassman- 
nian [SW85], etc.), study its consequences and better understand the underlying ge- 
ometry. 
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A Proof of Conjecture 7.1 up to o{l/N) 

Proposition A.l tlj{zi,Z2) is self-replicating at least up to o{l/N): 

^ 6^ij{zi, Z2) + ^(21, z)i){z, Z2) = o{l/N) 
Proof. Let us start from 11)12 = il^{zi, Z2) written as 



with 



1 Q12 I I f , ,(0) , 1 ®'l2 / /" f , ,(0) 



T^)Fi + UTr^-^)K" (A.2) 



-012 e/ ' 6V012 

where, to shorten notations, 1 means zi, 2 means Z2, § means the contour integral 
around ^B-cycles (indices are understood in tensor notations, i.e. contracted with the 
indices of derivatives of 0), and 0i2 means 

612 = 0(wo + 227r(u(zi) -u(z2)) |r), = 0(wo | r) (A.3) 

We remind that (Eqn. 5.3): 

(0). X B{zo,z)B{zuz)B{z2,z) 

= \ ^^l dx(.)dy(.) ^^-'^ 

and by special geometry: 

F- =m = y Res ^^^^^ (A.5) 

where we have denoted for short 

dv{z) = 2i7T du{z) = 1 B{z, ■) (A.6) 

Jb 
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The expression for 1^12 is thus: 



12 



1 rl rl 




(0) 



1 9i" 
26 



6' 



12 



(0) 



1/6' 



12 



6'\ 



612 6r^ 



(1) 



6' 



12 



6 



12 



2V6 
1/6^ 
6V612 

eJ 



12 



1 A-i 




(0) 



6^ 



^0 



+ Res 



dX{z)dY{z) 



1/6' 



- 



12 



+- 



1 /6'' 



12 



2V6 



12 



(dv(z))M5i2(^) 



1/6 



12 



2V612 
6"' 



6 



(dv(.))^ 



(A.7) 
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We need to apply the insertion operator jj 6^ to ?/'i2, and obtain the result up to o(l). 
Only the variation of NFq appearing in the Theta functions contributes to this order: 



TV 



12 



6' 



12 



6'\' 



612 6 
+ Res 



dX{z')dY{z') I2V6 



1/6; 



12 



112 



l\\2 



2 V 612 

+0(1) 

This allows us to compute: 
1 



1 /Pi'" 



(dv(^'))^d5l2(2') + 7 



dv(^')(d5i2(^')) 

" 6"' 



6 V6 



12 



6 



(dv(^'))n dv(2 



(A.8) 
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2iV 

2iv've 



+ ^ 5.^^12 + o(l/iV) 



(d5i2W)'+f|i^)dv(^')d5i2(^') 



iV ^z'^a, dX{z')dY{z') 1 2 
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i / KJ12 



^2V6i2 
d5i2(^) 4 
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6 



(dv(.'))^ + 



Ar2 



612/ 
5^2 + o(l/iV) 



612 6 



dw{z) 
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F[ dw{z) 



+— y Res 



^dX{z')dY{z') 1 2 
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dv(^')(d5i2(^'))'dv(z 



2 ve 
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12 







'(dv(/))'dv(z)} 



(A.9) 



We may transform the first line using the refined duality equation established in The- 
orem 3.3, which is a consequence of the Fay identity satisfied by the Theta function of 
the spectral curve. It can be rephrased as: 



d5i2(z) 



ei2 e 
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^1.^.2 eei2 
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On the other hand we have 
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Let us now compute: 
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1 /P)W p./// o/// p./// 1 

+-,i?r + ?f-?r- 7r)(dv(-'))^ - Usu^')dsM)ds,,{z') 

+o(l/iV) (A.12) 

The coefficient of F[ vanishes, as we can see by computing the gradient of Eqn. A. 10 
to leading order in N (recall that enters in the definition of our through the point 
Wo = ATF^): 

(0^ - ) dv(.) 

^12 61,62. /e;, , 6^2 6;^ 6' 



+ r£^_ri2_r^ ^A.i3) 



Ei,E,2 6612 V61, 6,2 612 6 
This identity can also be seen as a consequence of Fay identity. 
Let us now study the residue term of Eqn. A.12, which we write 

First, notice that by construction, Hi2{z',z) is a meromorphic 1-form in the variable 
z, which means that it has trivial monodromy when z goes around a non-trivial cycle. 
It may have simple poles at 2 = 2;i or 2; = 22 coming from the ratio of prime forms, 
but the expression in [■ ■ ■] vanish when 2 = 2:1 or Z2, so -^12(2, z') is actually regular at 
2 = 2:1 or ^ = ^2 • It may also have a singularity a.t z = z' coming from the term [■■■], 
which is atmost a double pole. To leading order when z ^ z', we find: 



+0{CMz)) (A.15) 

where is our notation for a local coordinate centered at z'. We use again Eqn. A. 10 
and find: 

HM,z) = idS,,(j'f +(|a)dvMdS,2M + i(||-|^)(dvM)' 

+ok;'(z)) 

+0{CMz)) (A.16) 
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We can rearrange the terms: 



0'l2 0' 



12 



e 



12 



e 



ei2 e 

dv(^' 
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and according to Eqn. A. 13, we see that Hi2{z',z) E 0{^~^{z)). Therefore, Hi2{z',z) 
is a meromorphic function whose only singularity is a pole atmost simple at z = z' . 
But a meromorphic function cannot have a single simple pole, so Hi2{z,z') must be 
holomorphic, and we can write it: 



H^2{z,z') = h2{z') dv(z) 



(A.17) 



Since the prefactor hi2{z') is independent of z, we may compute it by specializing to 
z = Zi'm Hi2{z,zi) defined from Eqn. A. 12. Doing so, we obtain: 
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All the terms eventually cancel each other. Thus Hi2{z,z') = 0, and coming back to 
Eqn. A. 12, this proves: 

— 6Jmlj{zi,Z2) + ^ = o{lN) A.19 
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